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THE MOTION OF A SATELLITE OF THE MOON 


by 

Giorgio E. O. Giacaglia. James P. Murohv, 
and Theodore L. Felsentreger 

in collaboration with 

E. Myles Standish, Jr,, and Carmelo E. Velez 


SUMMARY 


An analytical solution to the problem of the motion of 
a satellite of the moon is presented. Perturbations of short 
period and of intermediate period are derived through the 
application of the von Zeipel method. Long period pertur- 
bations are obtained through the use of elliptic integrals. 
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THE MOTION OF A SATELLITE OF THE MOON 

by 

Giorgio E. O. Giacaglia, James P. Murphy, 
and Theodore L, Felsentreger 

in collaboration with 

E. Myles Standish, Jr., and Carmelo E. Velez 
1. Introduction 


The motion of a satellite of the moon, or lunar orbiter, is analyzed. 
The solution is developed in powers of 10 -2 . Thus, a first order quan- 
tity is of order 10 -2 , a second order quantity is of the order 10~ 4 , and 
so on. The Hamiltonian for the "main problem" consists of zero, first, 
and second order quantities. The higher order Hamiltonian is of third 
order and consists of two parts: the first part contains terms generated 
by coupling of lower order terms and the second part consists of terms 
added by considering further perturbing forces, such as solar radiation 
pressure, physical libration, non- sphericity of the Earth's potential 
field, the attraction of the sun, etc. Additional terms are also produced 
by considering the eccentricity and inclination of the moon's orbit. 

In order to retain the relative orders of the disturbing forces, it is 
necessary to restrict the semi-major axis of the orbiter to about four 
moon radii or less. Further, certain restrictions must be made on the 
eccentricity and inclination in order not to invalidate the solution. 
Therefore, the following assumptions on the semi-major axis, eccen- 
tricity, and inclination are made: 

a < 4 moon radii 
.01 < e < .75 
| sin 1 1 > .01 . 

The small parameter of first order is n| , which is the mean motion 
of the moon's mean longitude . The small parameter s of second order are 
J 2 > J 22 > and (n c /n) 2 where J 2 and J 22 define the principal part of the 
oblateness of the moon and (n c /n) 2 is the square of the ratio of the mean 
motion of the moon to that of the orbiter. The small parameters of 
third order are J 3 , J 4 , J 5 , j 2 (n ( /n) 2 , (n # In) 2 , (r^/n) 2 sin(i c /2), (n e /n) 2 e 4 , 
07 » ( n « /n ) 3 , and an 0 /n. The quantities J 3 , J 4 , and J 5 are higher oblateness 
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parameters of the moon, and j 2 is the principal term in the oblateness 
of the earth. The quantity (n $ /n) 2 is the square of the ratio of the 
mean motion of the earth to that of the orbiter. The fact that the moon's 
orbital plane is inclined to its equator and the fact that the moon's orbit 
about the earth is elliptical give rise to the two small parameters of 
third order (n c /n) 2 sin (i c /2) and (n 5 /n) 2 e c , respectively. The radi- 
ation pressure gives rise to the third order parameter a. (n t /n) 3 is 
the cube of the ratio of the mean motion of the moon to that of the or- 
biter. Finally, an 0 /n is the correction due to physical libration. 

For the moon the values of J 2 and J 22 published by Jeffreys 
(Reference 5) are adopted, i. e. 


J 2 = 2.41 x i<r 4 
J 22 = 0.21 x IQ- 4 


The longest meridian of the moon contains the line joining the centers 
of mass of the earth and the moon. The right-handed, rotating, seleno- 
centric coordinate systems adopted for this problem will then be as 
follows: The z-axis is the rotational axis of the moon, and the xy -plane 

is the moon's equatorial plane. The x-axis passes through the moon's 
longest meridian, and is assumed to rotate with the motion n* . 

For a semimajor axis of 4 moon radii, 

2 4.8 x lO* 4 


(^) - 2.7 x lO' 4 

so that the oblateness terms in the lunar potential and the earth pertur- 
bations are both about of the same order. The largest oblateness term 
of the earth is then of the order of 10~ 6 . The perturbation of the sun is 
of third Qrder. 

Another perturbation to be considered is the effective radiation 
pressure of the sun, whose strength is about 1 X 10“ 4 dyne/cm 2 . If the 
area-mass ratio of the orbiter is 1.5 X 10 -1 cm 2 /g, then the disturbing 
acceleration due to radiation pressure is also of third order. 
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It is plain to see that a theory which includes the moon's oblateness 
and the earth perturbations should also include the solar perturbations 
in the "main problem" if the semi-major axis is above four moon radii. 

2. Equations of M otion 

The first step is to determine the equations of motion for the grav- 
itational fields of the moon as a primary and the earth and sun as per- 
turbations. 

The following notation will be used: 

subscript 0 : moon 

" 1 : orbiter 

" 2 : earth 

" 3 : sun. 

In an inertial system the equations of motion are 

m j p. = grad p . U (j = 0, 1, 2, 3) 

where 

,, ,2 / m 0 m i m 0 m 2 m 0 m 3 m l m 2 m l m 3 m 2 m 3\ 

\ r 01 r 02 r 03 r l2 r 13 r 23 / 

and p. is the radius vector of any one of the 4 bodies. If 77 ., £ . are 
the rectangular inertial coordinates of one of the bodies, then the equa- 
tions of motion can also be written 

*n (j = 0. !. 2, 3). 

Similar expressions hold for 7 ). and Cj. 

It is now convenient to refer the orbiter to a moon centered system, 
the moon to an earth centered system and the sun to a system whose 
origin is at the center of mass of the earth-moon system. Therefore 
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X 1 = - ^0 


x 0 = ^0 - ^2 


X 1 ~ ^ 3 ~ ~ ^3 


m 0^0 + m 2 ^2 
m Q + m 2 


Similar expressions hold for jr and z.. 


The equations of motion must be transformed accordingly. The partials 
in Eqs. (2) are then computed with respect to the new variables by mak- 
ing use of 


bu _ bu 3x j 

3 <f k Z_j3x c><f k 

* j=0 1 K 


(k= 0, 1, 2, 3) 


It follows that 

BU _ BU _ _BU m o BU_ 

B £ 0 Bx 0 Bx i m o + m 2 Bx 3 

BU _ BU 

Bxj 

BU _ BU m 2 BU 

B ^ 2 ~ Bx 0 m 0 + m 2 Bx 3 

_ _ 3 u_ 

B ^ 3 ~ Bx 3 


Since 

•• •• •; 

X 1 = - ^0 

it follows that 

x = (1 + Bu 1 1 BU 

x i \ m i + m 0 / Bx x m 0 Bx 0 + m Q + m 2 Bx 3 


The force function U must now be expressed in terms of the new co- 
ordinates. We have 

r oi = (£» - ^ 0 ) 2 + (Vi - v 0 ) 2 + (Cj - £ 0 ) 2 = x i + yi + z i = r i 
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r 02 = <fo - # 2 ) 2 + • • • = *3 t • • • = rj 


r 03 = (f> - f«) 2 + • • • 




x 3 + y 2 + z l + 


/ m 2 \2 2m 2 

< x o + y<> + z o> (x o x 3 + y 0 y 3 + z o z 3 ) 


2 ( m 2 V 2 2m 2 

r 3 + \ m 0 + m 2 y r ° m 0 + m 2 r o r 3 


r l2 = (^1 * ^2> 2 + * • * = ( X 1 + x 0> 2 + • • • = r i + r 0 + 2r l 


r L = <£i - ^ 3 ) 2 + 




*f . . . — 


= r\ + r 2 - 2r x . r 3 + f - 2 -- \ 
13 13 \ m o + m 2 / 


m„ \ 2 2m 2 

r o + i^Tl^ ( r o * r i - r o * ' 3 ) 


r 23 = <£2 - ^3> 2 + . • • 


= (- X 3-S^ X o) 


2 / m 0 \ 2 2 2m 0 
= r‘ + ( 1 r i + r„ • r 


3 \ m Q + m 2 ) 0 m 0 + m 2 0 * 3 


Since U contains the inverses of these radii, note that 


J_ _ J_ 

r oi ~ r i 


j_ _ 

r 02 r & 
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- 1/2 


r 03 


1 



2m 2 
m 0 *}■ 


L o 

— cos 
r. 



r l 2 


1 





COS 



1/2 



2m 2 

+ 

m 0 + rn 2 



cos s i0 


r cos S 03 

r 3 


1/2 



The angles used are shown in Figure 1. For obvious reasons the 
angle S Q3 is expressed in terms of Sq 3 , as follows: 


Since 


cos S' 


03 


r 0 • r 3 


r 0 r 3 


and 




r n + r~ = r, 

m 0 + m 2 0 3 3 ’ 


then 




rl + r n r, cos S 


cos 


c , „ m 0 +m 2° 03 
- — 


03 


03 


r 0 r 3 


m n 


r; m o + m 2 r ' 


+ — cos S 


03 
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/ 

/ 



or, finally 


Certainly, 


COS S °3 = r 3 ( cos S 03 
to third order, r 3 ~ r 3 . 


r o m 0 \ 

" m o + m 2 / ' 

so that 


cos S Q3 = cos 



m n 


m Q + m 


2 


On the other hand, S 10 may be replaced by 180° - S' 0 where Sj Q is 
the selenocentric elongation between the earth and the orbiter. 

3* Oblateness Terms 

If the plane of reference is the lunar equatorial plane, then the 
disturbing force per unit mass may be written as 
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u, 


OBL. 




n= 1 


J P (sin fi) cos m (X-A. ) 

nm nm y v nm 7 


The only terms to be included are J 22 , J 20 = - J 2 , J 30 J 3 » J 40 J 4 » 

andj 5 Q = -J 5 . Therefore, 


U 


OBL. 



J P 

J n n 


(sin (3) 


~ (^) J 22 P 22 < sin ft cos 2 ( V - X 2 2 )| ’ 

The angle /3 is the latitude of the orbiter with respect to the equator of 
the moon, A.' its longitude reckoned from any fixed direction and A 22 the 
longitude of the moon's longest meridian from the same fixed direction. 
Note that X! and /3 can be expressed in terms of the coordinates x x , yj, z\ 
of the orbiter. However, A .22 will contain the time explicitly, since 


k 22 = ^22 C°) + t 

where y is the frequency of rotation of the moon around its axis. 
Further, if we neglect the physical libration of the moon, the longest 
meridian is always pointing toward the earth and y t ~ n* . 


4. Canonical Equations and Gravitational Terms 

If we choose as canonical variables the Delaunay set 


L =y/i. Q a, -t = mean anomaly, 

G = LVl - e 2 , co = argument of pericenter, 
H=GcosI, 0 = longitude of ascending node. 


where a,e , and I are the semi-major axis, eccentricity, and inclination 
respectively, and where /i 0 = k 2 m ol where k is the Gaussian constant and 
m Q is the mass of the moon, then the equations of motion become 


L = 


BF 

3/ ’ 




/ 




ft = - 


o-. 

BF 
BH ’ 
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where 


'■v pi 

F = ^£2 + U GRAV. + U OBL. 


and where U GRAy . is a function which has to satisfy the equation 


Xl = Bxj ( U(3RAV - + r,°) 


The force function U depends on the new variables x through the relations 


r 2 — x 2 + v 2 4 - z 2 

i oi x i + y 1 + z i 


r 02 = X 3 + *0 + Z 0 


r 2 - 
L 03 


= ( x 3 - *°) 


r i 2 = < X 1 + x o)^ + 


: 13 = [( X l - X 3> + 


m„ 


-,2 


m 0 + m 2 o A 


23 


(- 


rl, = I - x, - 


m n 


3 m Q + m 2 0 


+ 


Then, 


BU , 2 J x o 

— - k 2 m Q m 2 — - m Q m 3 

r 02 


X, X n 

3 m 0 + m 2 0 


1 03 


(-—) 

\ m 0 + m 2 / 


m„ 


x i + x o 


X, - X, + 


- m i m 2 7 ~ 

r i 2 


~ m i m 3 


V 1 1 

1 3 m Q + m 2 0 / m, 


13 


(—) 

\ m 0 +m 2 / 
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Then, 


_3_ 

Bx, 


,2 r 3 ’ F 1 . 2 m 0 m 2 m 3 r 0 ‘ r l 

k m 0 m 3 1 + k 


1 03 


m ° + r o 3 


i 2 1 i 2 r 3 ' r i , 2 m 0 m 2 m 3 r 0 ‘ r l 

k m, m, k 2 m 2 m, k 2 

1 3 r i3 2 3 r 3 


23 


m n + r 3 

Z r 2 3 


:: m i + m 0 3 f t 2 m 0 m i 2 1 l2 1 1 

- - — < 4- m m, + k^ m. m, > 

1 m i m 0 3 Xl [ r 01 1 2 r 12 * - r 13 j 


1 B J 1 2 r O * r l , 2 m 0 m 2 m 3 r 3 ' r l 

^ — < - k 2 m 0 mj + k 2 


m o Bxj 


02 


m 0 + ^ r o 3 


m n m 0 m, • r 1 ^ 

+ k 2 m x m 2 

r 12 


j 2 “0 2 3 *0 M 


(m Q + m 2 ) 2 r 


03 


+ j^2 m l m 2 m 3 1 m ° m 2 m 3 r 3 ’ r i 

m 0 + m 2 r 13 m 0 + m 2 r 3 3 


2 m 0 m 2 m 3 r 0 * F 1 

- k 2 >• + 

( m 0 + m 2 ) 2 r 23 


1 B 


m o + m 2 B X j 


k 2 m 0 m 3 ^ + k 2 m ° m2 m3 F ° ‘ 

r 3 


r 3 (m 0 + m 2 ) r 3 


03 
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Therefore, 


where 


and where 


,2 1 i 2 r 3 • r i 

k 2 m. m- k z m 9 m, 

1 3 r i3 r 3 

13 r 23 


- k 2 


m 0 n *2 m 3 F 0 * r i 

m n + m 9 3 

U r 23 


-n I m n + m. m 0 

= -i. J k 2 - 1 + k 2 -1 + k 2 -L + k 2 m. 


1 r oi 


12 1 13 


r o ' r i 

2 r 3 

1 02 


i 2 //I 

- k m 3 r i • r 3 T 


1_ 

3 

' 03 


k 2 m k 2 m k 2 

U. . = + + + k 2 m. 

GRAV. r r r 2 3 


r o • r l 


01 12 ‘13 


02 


- k 2 


m 


F 1 • f 3 

3 r 3 
x 03 


m 9 

// _ ^ 

r , r q "■ r a 

3 3 m Q + m 2 0 


oi 1 1 


r 02 = f 0 


13 


03 


series in Legendre Polynomials 
series in Legendre Polynomials. 
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The term k 2 mj/r 01 can, of course, be neglected. Finally, 


F = 




^3 

+ — + M 2 

r l 3 


^3 F 1 


03 


+ 


^0 



(sin / 3 ) + 



(sin yd) cos 2 (A.'- k 22 


Now 


-r 0 r, cos S; Q 


r x 

= - — cos s; 0 


and 


Then 


1 1 r i 

— = r + ~ cos s io + 
r 1 2 r o rj 


P=2 


p (cos s; 0 > 



^2 


12 


+ 




(cos s; 0 ) 


and since M 2 /r 0 is independent of the position of the orbiter. 


F = 




2L 2 





P p ( cos SJo) + 


13 





(s in yS) + 



(sin y6) cos 2 (A.'-\ ) 


13 


A further simplification can be made if one considers the following 
expansion: 


r 


2 

13 





= T l + r 03 * 2 r i r 03 005 S 13 = r 03 


/ r i ^ F 1 

1 + f — ] -2 — - cos S" 3 

\ r 03/ 03 


l-= '-E (li-T p p s; 3 


13 ‘03 " \ 0 3 ; 


1 r i COS ^13 1 


1 03 


03 


— H (—) P P (cosS i'3' 


03 “ \ 03 


Also, 


r i ' r 3 r i r 03 005 S 1 3 r i COS S'; 


13 


03 


03 


1 03 


Since r Q3 is independent o£ X x , y x , z 1$ 


P = 2 


^3 y -1 /^i_v 

r 03 4^ '03/ 


p„ (cos s; 3 ) 


^0 



^ R -\2 


— J J n p n (sin + (~J J 22 P 22 (sin /3) cos 2 (V- 


K 


22 
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From now on the subscript 1 is omitted. 


5. Analysis of Terms 

Next the order of magnitude of the gravitational terms and the 
oblateness terms is computed. 

Consider the factor 


Then 


M 0 Mo 


r 3 a 3 (1 - e 2 ) 3 


+ periodic terms in f 


wnere 


n 2 a 3 = fi Q 


so that 




= n 2 , 




R 2 

R 2 = n 2 1 


Mo 




a 3 (1 - e 2 ) 


2n3 


(1 - e 2 ) 3 L 2 (1 - e 2 ) 2 


* US)']. 


2L 2 (1 - e 2 ) 


2n3 L va 


The relative magnitude of the greatest oblateness terms are then given 
by 


2 



J 2 


and 



J 


22 ’ 


Consider a very low orbiter for which a = 2R S . Then 
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2 (^“) 2 J 2 - 2x \ x 2x lO" 4 = 10- 4 
2(7) J 22 - 2 x \ x 0.2 x lO’ 4 = lO' 5 . 

Consider the earth's perturbation. It is dominated by the factor 

— (L 

r o \ r 0 

I 



or, neglecting periodic terms in f, 

„ 2 aUl-e 2 ) 2 ^( 1 -e 2 ) 2 C 

4 ‘ ^ I 




Since is small (~ 1/80), the relative order of magnitude of these 

terms is given by 2 (n 0 /n) 2 . For a high or biter (a - 5 R^) , 



2 


- 2 x 1 ° -3 • 


Consider the sun's perturbation. The distance r 03 is not well defined 
since the motion of the moon around the sun is far from a Keplerian 
motion. Instead of r Q3 , r 3 should be used. Then 
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The terms due to the sun's perturbation are therefore 
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(cos S" 3 ) 



m 2 \ 
m 0 + m 2 / 


cos S 


03 



* 


+ P 3 ( COS S 13^ 



COS S 03 + 



+ 


and they are dominated by the factor 


V - 3 a 2 (1 - e 2 ) 2 n 3 a 2 (1 - e 2 ) 2 /x 2 (1 - e 2 ) 2 

~ ^ ~ ' 7 0 + M 2 ^ 

3 1 + 

^3 

The relative order of magnitude is given by 

2 (?) " 2 * 10 ’ 5 

2 

for a high orbiter. Since, in the extreme case, a/a 0 - 3 X 10" , and if 
terms whose magnitude is less than 10 -5 are excluded, the disturbing 
function due to the earth is reduced to 



- (-i |p 2 (cos s; 0 ) +.L p 3 (cos s; 0 )l. 
r o \ r o/ ^ r 0 J 


For the sun, the only term to be included is 



P 2 (cos S" 3 ) 


since, in the extreme case, a/a 3 ~ 0.7 X 10“ 4 . 

Therefore, including terms of magnitude 10~ 5 , the Hamiltonian 
becomes 
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(^) : { P2 <cos s ‘ o) + ^ p 3 <cos s !« ) } 


77 (rj Pj (cos s “' ♦ t f £ £)’ * p n < si " n 

' L n=2 




J 22 P 22 ( sin Z 3 ') cos 2 (*•' " ^22) 


22 ' 


It is assumed that J 3 , J 4 , and J s are of the third order. At this time 
it is worthwhile to point out that the inclusion of the indicated spherical 
harmonics in the potential field of the moon in the manner described 
ahovp is an assumption, since this field is not very well known. It is 
possible that in the actual case one or more of J 3 , J 4 and J 5 may be 
greater than J 2 or J 22 . Our basic assumption is that the moon approx- 
imately satisfies the hydrostatic equilibrium conditions. 


6. The Angles S^andS^ 


In view of the form assigned to the oblateness terms, the plane of 
reference is the equator of the moon. Therefore, the next step is to 
express the angles S' Q and S|' 3 in terms of the orbital elements of the 
orbiter, the moon, the earth and the sun, with respect to that plane. 
The geometry is shown in Figure 2. 

The explicit form of S{ 0 and S" 3 requires the solution of two 
spherical quadrangles. This is now done. 
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7. The Cosine of Sj 0 


If x, y, z and x e , y e , z e are the rectangular coordinates of the 
orbiter and the earth, respectively, and r and r 0 their selenocentric 
distances, then 


cos s; 0 = 


r r 


Consider 


An =n _n 

© 


V = f + CO 
V © ~ ^ + 


Then, using vector notation, in the equatorial (moon) system indicated in 
Figure 3, 


cos w 




sin y_ cos 


sin v_ sm t 


r = r 


(cos co cos An - sin co sin An cos I) cos f \ 

“ (sin co cos An + cos co sin An cos I) sin f 

(cos co sin An + sin co cos An cos I) cos f 

- (sin co sin An - cos co cos An cos I) sin f 
sin co sin I cos f + cos co sin I sin f j 


This enables us to compute cos S' Q . 


8. The Cosine of S" 3 


In exactly the same way the cos S" 3 is obtained by substituting O 
for © in the formulas in section 7. 
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9 . The Main Problem 


The following approximations can be introduced: 

a. The earth's orbit around the moon (or vice versa) lies on the 
lunar equatorial plane. The error introduced by this approximation is 
proportional to the sine of half the inclination of the lunar orbit to its 
equator (~6°41'), or about 0.06. 

b. The orbit of the moon around the earth is circular and the motion 
uniform. The error is proportional to the moon's eccentricity or about 
0.055. 

c. The sun's perturbations are negligible. The relative magnitude 
for a moderately high satellite is about 0.05. 

d. The mean longitude of the earth, k s , is equal to ^- 2 2* 


If these approximations are introduced, the precision of the disturb- 
ing function is not higher than 10 -4 . It will be called the "disturbing 
function of the main problem, " and it is given by 


F 




r 2 P 2 (cos S 10 ) 




r 



{- J 2 P 2 (sin /3) + J 22 P 22 (sin cos 2 (k' - \ e )} , 


where 


e 



In this case the cos Sj 0 is much simplified. Since i c - 0 and v 0 + can 
be replaced by k 9 , it then follows 

cos S' 0 = cos v cos (Cl - k 9 ) - sin v sin (Cl - k 0 ) cos I = s. 
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The equations of motion are given by the canonical set 


BF 

• 

BF 

b7' 

l = 

" BL 

BF 


BF 

- — - ) 

00 = 

— 

dco 


BG 

BF 

n = 

BF 

BQ ’ 

~ BH 


Since the variables fi and A. m appear only in the combination fl - A. , where 
= n £ t + const, (approximation (b)), the degree of freedom is reduced by 
one by choosing as a new variable 

h = n - X . 

© 


The Hamiltonian must be modified accordingly, for 


h = Cl - n* = - — (F + n* H) 

c BH c 


(i.e., by using F + n * H in place of F). 


The Hamiltonian is still time dependent through \ 0 = n* t + const. 
Since the longest meridian is always pointing toward the earth, it is 
possible to choose the rotating system whose x-axis passes through this 
meridian. The final form of the Hamiltonian for the main problem is 
therefore 


n i r : 




Mo / R \ 

— {- J2 P 2 ( sin & + J22 P 22 ( sin & COS 2 ^} ’ 


where 


k = \' - K. 
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The equations of motion are (using in place of 


BF 
B/ ’ 


Bh 

BF 

~ BL * 

BF 

g 

BG 

h = -*I 

BH 


In the discussion that follows the solution of this system will be 
given using the already well known von Zeipel's method. 

10. Development of the Disturbing Function 

The notation used will be as shown in Figure 4. 



the earth 



and 


^ T 

cos Sj 0 = s = cos (v + fi - A. e ) - 2 sin 2 — sin v sin (\ e - ft) 


P 2 (s)^s 2 -4 = 


= y jcos 2 f cos (2g + 2h) + sin 2 (g + h) 


- sin f cos f sin(2g + 2h) + 4sin 4 ysin 2 h [cos 2 g - cos 2 f cos 2g 


+ sin f cos f sin 2 g ] + 4 sin 2 -isinh [sin f cos f cos (2g + h) 


- cos g sin (g + h) + cos 2 f sin (2g 


+ h)] }- 


1/2 . 


Further, 

sin/3 - sinl sinv = sin I sin (f + g) 

cos /3 cos (A. - h) = cos v 

cos /3 sin (\ - h) = cos I sin v , 

and therefore 

co s/3 co s\ - cos vcosh-cos I sinv sinh, 

from which 

P 2 (sin /3) = ysin 2 /3 -y = -|-sin 2 I sin 2 (f + g) -y 
P 22 (sin/3) cos 2A. = 6 cos 2 /3 cos 2 A. - 3 cos 2 /3 = 


6(£ 2 cos 2 f + x 2 sin 2 f + 2(y sin f cos f) - 3(1 - sin 2 I sin 


where 


* £ = cos g cos h - cos I sin g sin h 


X = - sin g cos h - cos I cos gsinh. 


Therefore, writing jjl c for /x Q , we have 


fit o n" r f 

F = + n * H + 1 cos 2 f cos (2g + 2h) 

2L 2 2 6 l 


+ sin 2 (g + h) - sin f cos f sin (2g + 2h) + 4 sin 4 sin 2 h [cos 2 g 


- cos 2 f cos 2g + sin f cos f sin 2g] + 4 sin 2 -^-sin h [sin f cos f 


• cos (2g + h) - cos g sin (g + h) + cos 2 f sin(2g + h)] j 

1 2 2 ^ f 3 T *2 t •2/r \ 

r + — r sin 1 sin (f + g) 

+ i J 2 + J 2 2 [ 6 ^ 2 cos2 f + X 2 sin 2 f + 2Cx sin f cos f) -3(1 - sin 2 1 sin 2 

= F q + Fj + F 2 , 
where 

M l 

F 0 = (Oth order) 

2L 2 

F x = n£ H ( 1 st order) 

F 2 = F - (F x + F 0 ) (2nd order) . 
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1 1 . Elimination of Short Periodic Perturbations 


The terms in F which depend on / correspond to the short periodic 
perturbations. According to von Zeipel's method, the elimination of 
these terms corresponds to the solution of the system 

F 0 '(L') = F 0 (L') 


Fj(H') =F 1 (H') + — i -2. 

dl BL' 


— j ' J 2 “'0 

** o + + “ — 

9h BH BL' 

ipVf ffo. 

+ 2 \ 3/ / 


where the generating function of the transformation 


(L, G, H, /, g, h) 


(L\ G', H\ /', g', h') 


is 


S = L'/ + G'g +H'h + Sj + S 2 + ••• , 


and the new Hamiltonian 


f' = f; + f; +f* + 


should be independent of / * . 
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In order to accomplish this last requirement, a particular solution 
is given by 


S x = 0 


F' = F 

r 2 2s 


!?» Hi. _ F 

3 / B L » 2 p * 


where F 2s and F 2p are, respectively, the parts of F 2 independent and 
dependent on 

If S 2 is determined from this system, then the relations between old 
and new variables are 


3S 2 

3S 2 

3S 2 

/' = / +-4, 

g = g + > 

h' = h +— 4, 

BL 

BG 

BH' 

BS, 

, 3S 2 

, 3S 2 

L = L' + 2 , 

G = G + 2 , 

H = H + — — 

3/ 

^g 

Bh 


Next, the "secular" part of F 2 is determined. By definition, 


t 2v 

F ”=5Tl F ’ d '- 


i-f 

277 I 

•'A 
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The integrals needed are 


±r 

2 77 J 

*'a 


r 2 cos 2 f d/ = (1 + 4 e 2 ) — 

^ 2 


± r 

27 , J 


r 2 d/ = ( 1 + — e 2 ) a 2 


r 27T 

^ f rJsi 


in f cos f d/ = 0 , 


which are evaluated using the relations (E = eccentric anomaly) 

d / = (1 - e cos E) dE 
r cos f = a (cos E - e) 

r sin f = a /l - e 2 sin E 
rra(l-ecosE). 


Furthermore, the following integrals are also needed: 


jl r j_ 

277 X r 3 


cos 2f d/ =0 
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sin f cos f d/ = 0. 


£ 7 

J o r 


These are evaluated using the relations 


dl =■ 


a 2 (1 - e 2 ) 1/2 


df 


1 1 + e cos f 


r a(l-e 2 ) 


If the preceding results are used, it follows that 


F' 


£ . -*u' 


2L 


>2 


+ n*H' + 


n * L ' 4 l/ 5 -3 G ' J 


16 6 


'2 , 


■)[(■* -S) 


3 (‘-^) ”■ 2h '] + 15 (‘ -£) [i (‘ ♦{?)’ 


cos 2 (g' + h') 


+ l 1 “ COS 2g< + T V 1 " G 7 ) COS 2 (g> " h ' } 


’]} 




and since 


3F' 

3/' 


0, 


then 


L* = const. 


The generating function and the short period perturbations are ob- 
tained next. Consider F in the form 


R“ 

F = 2I^ +n ® H+: ^7 i [- J 2 P 2 (S^/ 3 ) +J 22 P 22 ( sin P) COS 2 k] 


n: 

(t £ 


+ — i feF v>- 


and write 


U x = r 2 P 2 (s) 


2tt 

"»■-£ 1 

•'n 


r 2 P 2 (s) d/ 


U 


Ip 


r 2 P 2 (s)-U ls 


U 2 = _ P 22 (sin/3) cos 2\ 
r 3 


f 27f 


dl 
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The evaluation of the integrals is carried out using E or f as independent 
variables. 

It is found that 


K 


d l = - D ( — sin 2 I 


, 1 \ 3 

" 2 j + 4 


K sin 2 I' 


(1 - e 2 ) 


* 2-V-3/2 


,'3 


/ 3 • 
— SI 

U 


sin 2 I' 



/ 


U 2p d/ = (f>x' 2 -3 + 3 sin 2 I' cos 2 g ) D 


+ [ 6 (£' 2 ~ x ' 2 )- 3 si n 2 I' cos 2g] Q 


+ [uC'x +3sin 2 I' sin 2,]J 
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2a' 3 


(1 - e* 2 ) _3/2 sin 2 I' cos 2h 


j o' 



cos (2g + 2h) - 6 sin 4 


1 ■ 2 i 

— sin x h cos 
2 


2g 


+ 6 sin 2 


— sin h sin 
2 



f 3 j 1 

+ B J— sin 2 (g + h) + 6 sin 4 — sin 2 h cos 2 g 
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-6 sin — 
2 


sin h cos g sin (g + h) - 


2 


) 


+ C sin (2g + 2h) + 6 sin 4 ^ sin 2 h sin 2g 


+ 6 sin 2 sin h cos (2g + h) j- 


|(2 + 3e' 2 ) [(-1 + 3 cos 2 I') +3 sin 2 I 


»* 2 ri 

_2 


+ 15 e' 2 — (1 + cos I') 2 cos (2g + 2h) + sin 2 I 


2 T' 



(1 - cos I') 2 cos (2g - 2h) 


} 


cos 2 


cos 2g 


where 



C' = cos g cos h - cos I' sin g sin h 
x' = - sin g cos h - cos I* cosgsinh 



cos 2 f' d/ 


»2 




sin E' 


+ 



sin 2 E' 


— sin 3E* 

12 



3e' 



sin E' 


+ — e' 2 sin 2 E' - - — sin 3 E 1 
4 12 


C = 


J 


.'2 


sin f' cos f' d/ 




cos E' - — (1 + e' 2 ) cos 2 E' 


e _ _ f 

+ — cos 3 E 

12 



d/ = a'” 3 (1 - e' 2 ) 3/2 [f' + e* sin f*] 
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r cos 2 f* 

J r' 3 


- o'-3 


dl = a 


(1 


_^- 3/2 r - f' + — sin f' 
|_ 2 4 


+ — sin 2 f' + — sin 3 f' 
4 12 


= f sin f' cos f' d = a * -3 (1 _ ^ 2 r 3/2 L cos f _ 1 cos 

Jr' 3 [44 


2 f' 


cos 3 f 

12 


I 



cos 2 ( f ' + g) 


dl = a 


'- 3 (I - e ' 2 )- 3/2 


1 sin ( 2 g + 2 f') 


and 


+ -77 sin ( 2 g + f') + 7 - s in (2g + 3 f ' ) 


2 ' ” ’ 6 

E' - e' sin E' = / 


(def ines E') 


r' = a' (1 - e' cos E') (defines r') 


J_ _ 1 + e' cos f' 
r' a' (1 - e' 2 ) 


(de fines f ' ) 


In the integrals above, any constant with respect to the variable of 
integration ( E' or f' ) has been neglected. At this stage, S 2 is com- 
pletely defined. After some simplification, the final result is 
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00 I H-k 


fjL R 2 j ( 

S 2 = — |-2 (1 - 3 cos 2 I') [e' sinf' + (f'-/)] + sin 2 l' [3 e' sin(f' + 2g) 

8 n'a' 3 (l-e' 2 ) 3/2 

+ 3 sin (2f* + 2g) + e' sin (3f' + 2g)] j 

A 6 J i 

^ — — \6 s in 2 I* [2(f' - / ) cos 2h + e' sin(f ' - 2h) + e' sin(f ' + 2h)l 

n'a' 3 (l-e' 2 ) 3/2 

+ (l-cosI') 2 [3e'sin(f' + 2g- 2h) + 3sin(2f' + 2g - 2h) + e' sin (3f* + 2g - 2h)] 

+ (1 + cos I*) 2 [3 e' sin (f' + 2g + 2h) + 3sin(2f' + 2g + 2h) + e' sin (3f' + 2g + 2h)]j 
n 2 q # 2 

+ — - — — — — fl2 [-2 (1 - 3 cos 2 I') (2 + 3 e' 2 ) + 6 sin 2 I' (2 + 3e' 2 )cos 2h 
384 en' *" 

+ 30 e' 2 sin 2 I* cos 2g 

+ 15e' 2 (l-cosI') 2 cos(2g - 2h) + 15e' 2 (l + cosl') 2 cos(2g + 2h)] • (E* - / ) 

+ 4 (1 - 3 cos 2 I') [9 e' (4 + e' 2 ) sinE' - 9 e' 2 sin 2E' + e' 3 sin 3E'] 

+ 6 sin 2 I' [-9 e' (4 + e' 2 ) sin(E' - 2h) +9e' 2 sin(2E' -2h)-e' 3 sin(3E' - 2h) 

- 9 e' (4 + e' 2 ) sin (E' + 2h) + 9 e' 2 sin (2E' + 2h) - e' 3 sin (3E' + 2h) 

- 15 e' {(2 + e ,2 )-2(l-e' 2 ) 1/2 >sin(E' -2g) 

+ 3 {(2 + e' 2 ) -2(l-e' 2 ) 1/2 (1 +e' 2 )} sin(2E' -2g) 

-e' {(2-e' 2 ) -2 (1 -e' 2 ) 172 } sin(3E' -2g) 

- 15 e' {(2+e' 2 ) + 2(1- e ' 2 ) 172 } sin(E' + 2g) 

+ 3 {(2 + e' 2 ) +2(1- e' 2 ) 172 (1 + e' 2 )} sin (2E' + 2g) 

- e' { (2- e' 2 ) +2(l-e' 2 ) 172 > sin(3E' +2g)] 

+ 3 (1 - cos I') 2 [-1 5 e' {(2 + e' 2 )-2(l - e' 2 ) 172 } sin (E' -2g + 2h) 

+ 3{(2 + e' 2 )-2(l-e' 2 ) 172 (1 + e' 2 )}sin(2E' - 2g + 2h) 

-e' {(2-e' 2 ) -2(l-e' 2 ) 1/2 } sin(3E' -2g + 2h) 
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- 15 e' {(2 + e' 2 ) + 2(1 - e' 2 ) 1/2 } sin (E' + 2g - 2h) 

+ 3{(2 + e' 2 ) + 2(1 - e' 2 ) 1/2 (1 + e' 2 » sin (2E' + 2g - 2h) 
-e' {(2-e' 2 ) +2(1 -e' 2 ) 1/2 } sin(3E' +2g-2h)] 

+ 3(1 + cos I') 2 [- 15 e' {(2 + e' 2 )-2(l-e' 2 ) 1/2 } sin (E* - 2g - 2h) 

+ 3 {(2 +e' 2 ) -2(1 -e' 2 ) 1/2 (1 +e' 2 )} sin(2E' -2g-2h) 

-e' {(2-e' 2 )-2(l-e' 2 ) 1/2 } sin(3E' -2g-2h) 

- 15 e' {(2 + e' 2 ) + 2(1 - e' 2 ) 1/2 } sin(E' + 2g + 2h) 

+ 3 {(2 + e' 2 ) + 2(1 - e' 2 ) 1/2 (1 + e' 2 )> sin(2E' +2g + 2h) 

-e' {(2-e' 2 ) +2(l-e' 2 ) 1/2 } sin(3E' +2g + 2h)]j, 


so that, for the short periodic terms, 

S = V l + G'g + H'h + S 2 

to second order. 

The short period perturbations are given by 



g 


BS _ 3S 2 

BG' 6 + BG' 



h + 


BS, 

bh' 
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G 


BS 

3g 



H = — = H' 

Bh 



It is easier to compute the partial derivatives with respect to the 
Keplerian elements a', e', I' and then compute them with respect to 
L' , G' , H' as follows: 

!!?-!!? 3a' ^ Be^ _ 0 L' ^2 , _G^_ ^2 
BL' Ba' BL' Be' BL' “V* Ba' ‘ e 'L' 3 Be' 


^S_ 2 _^S_ 2 Be' + ^B_r G' *S 2 H' 

BG' Be' BG' + Bl' BG' e 'L' 2 Be' + G' 2 sinl' Bl' 


BS 2 3s 2 BI* 1 3S 2 

BH' ~ BI' BH' ~ " G' sin I' Bl' ' 

Furthermore, it is important to note that 


BF 

Be' 



f' 


BE' 

Be* 



n E' 


Bn' _ 3 n^ 

Ba' " ~ 2 a ' 
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The partial derivative with respect to l does not need to be com- 
puted^ since 




®<c J 2 


U 3P+' 


J22 


= J_(F -F'). 
n' 


To simplify the formulae, the following notations are now introduced: 
B 2>1 - e' sin f' + (f' - l) 

B 22 =3e' sin(f'+2g)+3sin(2f' + 2g)+e'sin(3f' + 2g) 


B 22 j = 2(f' -/) cos 2h + e' sin (f' - 2h) + e' sin (f* + 2h) 


B 22 2 = 3 e' sin(f' + 2g -2h) +3 sin (2f' + 2g - 2h) + e' sin(3f' + 2g- 2h) 
B 22 3 = 3 e' sin(f' + 2g + 2h) + 3 sin (2 f ' + 2g + 2h) + e' sin (3f ' + 2g + 2h) 


B c j = 9e' (4 + e' 2 ) sinE' -9e' 2 sin2E' + e' 3 sin3E' 


B c. 2 = " 9e ' (4 + e' 2 )sin(E'-2h) + 9e' 2 s in(2E' - 2h) - e' 3 sin(3E'-2h) 

- 9 e' (4 + e' 2 ) sin (E' + 2h) + 9 e' 2 sin (2E' + 2h) - e' 3 sin(3E' + 2h) 

- 15e'{(2 + e' 2 )-2(l-e' 2 ) 1/2 } sin(E' - 2g) + 3{(2 + e' 2 )- 2(l-e' 2 ) 1/2 (1 + e' 2 )>. 
. sin (2E'- 2g) - e' {(2-e' 2 ) - 2(1 -e' 2 ) 1/2 } sin(3E'-2g)- 15 e' {(2 + e' 2 ) 

+ 2 (1-e' 2 ) 1/2 } sin(E'+2g) + 3 {(2 + e' 2 ) + 2 (l-e' 2 ) 1/2 (l +e' 2 ) }sin(2E' +2g) 

- e' {(2-e' 2 ) + 2(1 - e' 2 ) 1/2 } sin (3E' +2g) 
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B c 3 = - 15 e' {(2 + e' 2 ) - 2(1 -e' 2 ) 1/2 } sin(E'- 2g + 2h) +3{(2 + e' 2 ) 

-2(l-e' 2 ) 1/2 (l +e' 2 )} sin(2E'-2g + 2h) -e' {(2 - e' 2 ) - 2(1 - e' 2 ) 1/2 } • 
. sin(3E'-2g + 2h) - 15e' {(2 + e' 2 ) + 2(1 - e' 2 ) 1/2 > sin(E' + 2g - 2h) 

+ 3 {(2 + e' 2 ) + 2(1 - e' 2 ) 1/2 (1+e' 2 )} s in (2E' + 2g - 2h) - e' {(2-e' 2 ) 

+ 2(1 - e' 2 ) 1/2 } sin(3E' + 2g- 2h) 


4 = - 15 e' {(2 + e' 2 ) - 2 (1 - e' 2 ) 1/2 } sin (E' - 2g - 2h) + 3 {(2 + e' 2 ) 

- 2(1 - e' 2 ) 1/2 (1 + e' 2 )} sin (2E' - 2g - 2h) - e' {(2 - e' 2 ) - 2(1 - e' 2 ) 1/2 } • 
• sin(3E' - 2g- 2h) -15 e' {(2 + e' 2 ) + 2(l-e' 2 ) 1/2 } sin(E' + 2g + 2h) 

+ 3{(2 + e' 2 ) + 2(l-e' 2 ) 1/2 (1 + e' 2 )} sin(2E' +2g + 2h) 

-e' {(2-e' 2 ) + 2(l - e' 2 ) 1/2 } sin(3E' +2g + 2h) . 


With this notation, it follows that 


*S 2 


3a' 


n — (-2 (1-3 cos 2 I')B 2 1 + (sin 2 I') B, } 

l6 a'(l-e' 2 ) 3/2 


3 n'R 2 J 


22 


16 a'(l-e' 2 ) 2/2 


{6(sin 2 I')B 22>1 t (1-cos I') 2 B 22 2 t (1 + cos I') 2 B 22 3 > 


7 a> 

— — (12 [-2(1-3 cos 2 I') (2 + 3 e' 2 ) +6sin 2 I' (2 + 3 e' 2 ) cos 2h 


+ 30 e f 2 s in 2 I' cos2g + 15e #2 (l-cosl') 2 cos(2g- 2h) + 15e' 2 (l+cosI # ) 2 # 
.cos(2g + 2h)]-(E'-/)+4(l-3cos 2 l') B c>1 +6 (sin 2 I') B e>2 
+ 3(1 - cos I') 2 B £ 3 +3(1 +cosl') 2 B (E 
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H I ^ + + + • ro 100 




(1 -e' 2 ) 3/2 


sin I' cos I' {-6 B 21 + B 2 2 ) 


1 n ' R ®J 2 2 


4 ( 1 - e * 2 ) 3/2 


sin I ' {6 (cos I' ) B 22 j + ( 1 - cos I' ) B 22 2 -(l + cos I ' ) B 22 _ 3 ) 


1 n'a' 2 /n \ 2 

— — f -f sin I' {12 [- 2 cos I' (2 + 3e' 2 ) + 2 cos I' (2 + 3e' 2 ) cos 2h 


1 0 e' 2 cos I' cos 2g + 5e' 2 (l- cos I ') cos (2g - 2h) - 5e' 2 (1 + cos I') cos(2g + 2h)] 


(E '-/) +4(cosr)B cl + 2 (cos I') B £>2 + ( 1 - cos I ' ) B ffi>3 - (1 + cos I * ) B £ 4 > 


e'n' R|J 2 


(1 -e' 2 ) 


2x5/2 


{-2(1-3 cos 2 I ' ) B 2 j + ( sin 2 I ' ) B 2 2 } 


1 "'R 2 J 2 


+ — 


8 (l-e' 2 ) 3/2 


{-2(1-3 cos 2 1 ' ) sin f + sin 2 1 9 [3 sin ( f '+ 2 g) + sin (3f # + 2g)] } 


+ - 


1 n#R ^J 2 


8 ( 1 - e' 2 ) 3/2 


{-2(1-3 cos 2 I') [e cos f ' + 1] +3 sin 2 !' [e' cos(f ' + 2g) 


+ 2 cos ( 2f ' + 2 g) + e ' cos ( 3f ' + 2 g)] } f— — \ sin f 


r' G ' 2 


3 e '"' R e J 22 


(1-e' 2 ) 


' 2 \5/2 


{6(sin 2 I')B 22>1 + (1 - cos I') 2 B 22>2 + (1 +cosl') 2 B 22>3 } 



i n ' K J 22 

+ T (6 sin 2 I' [sin (f'-2h) + sin (f'+ 2h)] +(l-cosl') 2 [3 sin (f '+ 2g - 2h) 

8 (l-e' 2 ) 3/2 


+ sin (3f '+ 2g - 2h)] + (1 + cos I') 2 [3 sin (f'+ 2g + 2h) + sin (3f '+ 2g + 2h)] } 


3 n ' R lJ 22 

+ — {2 sin 2 I' [2 cos 2h + e' cos ( f 2h) + e' cos (f ' + 2h)] 

8 (l-e' 2 ) 3/2 

+ (1 - cos I') 2 [e' cos(f ' + 2g - 2h) + 2 oos(2f '+ 2g - 2h) + e' cos (3f ' + 2g - 2h)] 


+ (1 + cos I') 2 [e' cos(f'+ 2g + 2h) + 2 oos(2f ' + 2g + 2h)+ e' cos(3f '+ 2g + 2h)]} • 


V* o' : 


«' L ' 2 \ . .. 

„ sin f 

12 J 


1 n 


>2 /n„\ 2 


a~ 


— * j |24e' [-2(1-3 cos 2 I* ) + 6 sin 2 I ' oos 2h + 10 sin 2 I* cos 2g 


128 e 

+ 5(1- cos I ') 2 cos (2g - 2h) + 5(1 + cos I') 2 cos (2g + 2h)] • (E '- / ) 

+ 4(1-3 cos 2 I') [3(4 + 3 e' 2 ) sinE' - 6 e' sin 2E' + e' 2 sin 3E'] 

+ 2 sin 2 I' [- 9 (4 + 3 e' 2 ) sin (E'- 2h) + 18 e' sin(2E'- 2h) - 3e' 2 sin (3E'- 2h) 
-9 (4 + 3 e' 2 ) sin (E'+ 2h) + 18 e' sin (2E'+ 2h) - 3 e' 2 sin (3E'+ 2h) 

- 15 {(2 + 3 e' 2 ) - 2(1 - e' 2 ) 1/2 + 2e' 2 (1 -e' 2 )" 1/2 } sin (E' - 2g) 


+ 6 e' (1 - 2 (1 - e' 2 ) 1/2 + (1 - e' 2 )- 1/2 (1 + e' 2 )} sin (2E* - 2g) 

- {(2 - 3 e' 2 ) - 2(1 - e' 2 ) 1/2 + 2e' 2 (1 - e' 2 )- 1/2 } sin (3E‘- 2g) 

- 15 {(2 + 3e' 2 ) + 2(1- e' 2 ) 1/2 - 2e' 2 (1 - e' 2 )- 1 ' 2 } sin (E'+ 2g) 
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+ 6 e' {1 + 2(l-e' 2 ) 1/2 -(l-e ,2 )- 1/2 (1 + e' 2 )} sin(2E'+2g) 

- {(2 - 3e' 2 ) + 2(1 - e' 2 ) 172 - 2e' 2 (1 - e' 2 )* 172 } sin (3E* + 2g)] 

+ (1 - cos I') 2 [-15 {(2 + 3e' 2 ) - 2(1 - e' 2 ) 172 + 2e' 2 (1 - e' 2 )" 172 } sin(E'- 2g + 2h) 

+ 6e' {l-2(l-e' 2 ) 172 + (l-e' 2 )~ 1/2 (1+e' 2 )} sin (2E'- 2g + 2h) 

- {(2 - 3e' 2 ) - 2(1 - e" 2 ) 172 + 2e' 2 (1 - e' 2 )" 172 } sin (3E'- 2g + 2h) 

- 15 {(2 + 3e' 2 ) +2(1 -e' 2 ) 172 -2 e' 2 (1-e' 2 )" 172 } sin (E'+ 2g-2h) 

+ 6 e' U +2(1 -e' 2 ) 172 - (1 - e' 2 )* 172 (1 + e' 2 )} sin (2E'+ 2g - 2h) 

-{(2-3e' 2 ) + 2(1 - e' 2 ) 172 - 2e' 2 (1-e' 2 )- 172 } sin (3E'+ 2g - 2h>] 

+ (1 + cos I') 2 [- 15 {(2 + 3e' 2 ) - 2(1 - e' 2 ) 172 + 2e' 2 (1 - e' 2 )" 172 } sin(E'- 2g - 2h) 
+ 6e' (1-2(1 -e' 2 ) 172 + (1 - e' 2 )" 172 (1 + e' 2 )} sin (2E'- 2g - 2h) 

- ((2 - 3e' 2 ) - 2(1 - e' 2 ) 172 + 2e' 2 (1 - e' 2 )" 172 } sin(3E'- 2g - 2h) 

- 15 {(2 + 3e' 2 ) + 2(1 - e' 2 ) 172 - 2e' 2 (1-e' 2 )" 172 } sin(E'+ 2g + 2h) 

+ 6e'{l + 2(l-e' 2 ) 1/2 - (1-e' 2 )- 1 / 2 (1 + e' 2 )} sin(2E'+ 2g + 2h) 

- ((2 -3e' 2 ) + 2(1 - e' 2 ) 172 - 2e' 2 (1 - e' 2 )" 172 } sin(3E'+ 2g + 2h)]| 


+ (l 2 [-2(1-3 cos 2 l')(2 + 3e' 2 ) + 6 sin 2 l|* (2 + 3e' 2 ) cos 2h 

384 £ J ^ 

+ 30 e' 2 sin 2 I' cos 2g + 15e' 2 (l - cos I' ) 2 cos (2g- 2h) + 15e' 2 (1 + oosl') 2 cos(2g + 2h)] 


+ 12(1-3 cos 2 I ') [3 e' (4 + e' 2 ) cos E' -6e' 2 cos 2E' +e' 3 cos3E'] 
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+ 18 sin 2 I' [- 3e' (4 + e' 2 ) cos (E' - 2h) + 6e' 2 cos (2E'- 2h) - e' 3 cos(3E'- 2h) 

- 3e' (4 + e' 2 ) cos (E' + 2h) + 6e' 2 cos(2E'+ 2h) - e' 3 cos(3E'+ 2h) 

- 5 e' {(2 + e' 2 ) - 2(1 - e' 2 ) 1/2 } cos (E'~ 2g) 

+ 2 {(2 + e' 2 ) - 2(1 - e' 2 ) 1/2 (1 + e' 2 )} cos (2E'- 2g) 

-e' {(2-e' 2 ) - 2(1- e' 2 ) 1/2 } cos (3E'-2g) 

-5e'{(2 + e' 2 )+2(l-e' 2 ) 1/2 }cos(E' + 2g) 

+ 2 {(2 + e' 2 ) + 2(1 - e ,2 ) 1/2 (1 + e' 2 )} cos(2E' + 2g) 

-e' {(2-e' 2 ) +2(1 -e' 2 ) 1/2 > cos(3E'+2g)] 

+ 9(1 _ oos I') 2 [- 5 e' {(2 + e' 2 ) - 2(1 - e' 2 ) 1/2 } cos (E*- 2g + 2h) 

+ 2 {(2 + e' 2 ) - 2(1 - e ,2 ) 1/2 (1 + e' 2 )} cos (2E'~ 2g + 2h) 

- e' {(2 - e' 2 ) - 2(1 - e' 2 ) 1/2 } cos (3E'~ 2g + 2h) 

-5e'{(2 + e' 2 )+2(l-e ,2 ) 1/2 } cos ( E' + 2g - 2h) 

+ 2 {(2 + e' 2 ) + 2 (1 - e' 2 ) 1/2 (1 + e' 2 )} cos (2E'+ 2g - 2h) 

-e' {(2-e' 2 ) + 2(1 - e' 2 ) 1/2 } cos (3E' +2g - 2h)] 

+ 9 (1 + cos I' ) 2 [- 5 e' {(2 + e' 2 ) - 2(1 - e' 2 ) 1/2 } cos (E'- 2g - 2h) 

+ 2 {(2 + e' 2 ) - 2(1 - e' 2 ) 1/2 (1 +e' 2 )} cos (2E'- 2g - 2h) 

- e' {(2 - e' 2 ) - 2(1 - e' 2 ) 1/2 } cos (3E'- 2g - 2h) 
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- 5 e' {(2 + e' 2 ) + 2(1 - e' 2 ) 1/2 } cos (E'+ 2g + 2h) 

+ 2 {(2 + e' 2 ) + 2(1 - e ,2 ) 1/2 (1 + e' 2 )} cos (2E'+ 2g + 2h) 

- e' {(2 - e' 2 )+ 2(1 - e' 2 ) 1/2 }cos (3E'+ 2g + 2h)]} -ll sinE* 


3S n ' R 2 T 

— - = — 2 sin 2 I' (3 e' cos(f ' + 2g) + 3 cos (2f '+ 2g) + e' cos (3f'+ 2g)} 


3g 


4 (l-e'*) 3/I 


1 11 J22 

+ *(l-e' 2 ) 3/2 


{(1 - 00 s I') 2 [3e' cos (f' + 2g - 2h) + 3 oos(2f'+ 2g - 2h) 


+ e' cos(3f'+ 2g - 2h)] 

+ (1 + cos I ') 2 [3e* oos(f'+ 2g + 2h) + 3 cos(2f'+ 2g + 2h) + e' cos(3f'+ 2g + 2h)] > 
+ {” 60 e' 2 [2 sin 2 I * sin 2g + (1 - cos I ') 2 sin(2g- 2h) 

+ (1 + cos I') 2 sin(2g + 2h)] • (E'~ / ) + 2 sin 2 I '[15 e'{(2 + e' 2 ) 

-2(1 -e' 2 ) 1/2 } cos (E'-2g)-3{(2 + e' 2 ) -2(1 -e' 2 ) 1/2 (1 + e' 2 )} cos(2E'-2g) 


+ e'{(2-e' 2 )-2(l-e' 2 ) 1/2 } cos (3E'-2g) 

- 15 e' {(2 + e' 2 ) + 2(1 - e' 2 ) 1/2 } cos (E'+2g) 

+ 3 {(2 + e' 2 ) + 2(l-e' 2 ) 1/2 (1 + e' 2 )} 00 s (2E'+2g) 

- e* - e' 2 ) + 2(1 - e' 2 ) 1/2 } cos(3E'+ 2g)] 

+ (1 - cos I') 2 [l5e'{(2 + e' 2 )-2(l - e' 2 ) 1/2 } cos(E'- 2g + 2h) 

- 3 {(2 + e' 2 ) - 2(1 - e ' 2 ) 1/2 (1 + e' 2 )} cos (2E*- 2g + 2h) 
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+ e' {(2 - e' 2 ) - 2(1 - e' 2 ) 1/2 } cos (3E'~ 2g + 2h) 

- 15e' {(2 + e' 2 ) + 2(1 - e' 2 ) 1/2 } cos (E' + 2g - 2h) 

+ 3 {(2 + e' 2 ) + 2(1 -e' 2 ) 1/2 (1 +e' 2 )} cos (2E'+ 2g - 2h) 

-e' {(2 - e* 2 ) + 2(1 - e ,2 ) 1/2 } oos(3E‘+ 2g - 2h)] 

+ (1 + cos I 9 ) 2 [15 e' {(2 + e' 2 ) - 2(1 - e' 2 ) 1/2 } cos(E'- 2g - 2h) 

- 3 {(2 + e' 2 ) - 2(1 - e' 2 ) 1/2 (1 +e' 2 )} cos (2E'- 2g- 2h) 

+ e' {(2 - e' 2 ) - 2(1 - e' 2 ) 1/2 } cos(3E'- 2g - 2h) 

- 15 e' {(2 + e' 2 ) + 2(1 - e' 2 ) 1/2 } cos(E'+ 2g + 2h) 

+ 3{(2 + e' 2 ) + 2(l-e' 2 ) 1/2 (l+e' 2 )} cos (2E'+ 2g + 2h) 

-e' {(2-e' 2 ) + 2(l-e' 2 ) 1/2 > cos (3E'+ 2g + 2h)] j 


BS 2 1 n'RcJ 


C J 22 


3h 4 (1 - e' 2 ) 3 ^ 2 


{6 sin 2 I' [-2(f / ) sin 2h - e' cos(f 2h) + e # oos(f # + 2h)] 


- (1 - cos I') 2 [3 e' cos(f ' + 2g - 2h) + 3 cos(2f' + 2g - 2h) + e' cos(3f' + 2g - 2h)] 


+ (1 + cos I 9 ) 2 [3e' cos (f ' + 2g + 2h) + 3 cos (2f # + 2g + 2h) + e cos(3f '+ 2g + 2h)] } 


+ r^- n (— ~\ |36 [-2 sin 2 I' (2 + 3 e' 2 ) sin 2h + 5e' 2 (1 - cos I') 2 sin(2g - 2h) 


192 e 


- 5 e* 2 (1 + cos I') 2 sin(2g+2h)] •(£*-/) 
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+ 6 sin 2 I' [9 e' (4 + e' 2 ) cos (E'- 2h) - 9e' 2 cos (2E'- 2h) + e' 3 cos (3E'- 2h) 

- 9 e' (4 + e' 2 ) 00 s (E* + 2h) +9e' 2 cos (2E' + 2h) - e' 3 cos(3E'+2h)] 

+ 3(1 - cos I') 2 [- 15 e' {(2 + e' 2 ) - 2(1 - e' 2 ) 1/2 } cos(E'- 2g + 2h) 

+ 3{(2 + e' 2 )-2(l-e' 2 ) 1/2 (1 + e' 2 )} oos(2E'- 2g + 2h) 

-e' {(2-e' 2 ) -2(1 -e' 2 ) 1/2 } cos (3E'- 2g + 2h) 

+ 15 e' {(2 + e' 2 ) + 2(1 - e' 2 ) 1/2 } cos(E'+ 2g - 2h) 

- 3 {(2 + e' 2 ) + 2(1 - e' 2 ) 1/2 (1 + e' 2 )} cos (2E'+ 2g - 2h) 

+ e' {(2-e' 2 ) +2(l-e' 2 ) 1/2 } cos (3E'+2g-2h)] 

+ 3(1 + cos I') 2 [l5 e' {(2 + e' 2 ) - 2(1 - e' 2 ) 1/2 } cos(E'- 2g - 2h) 

-3{(2 + e' 2 )-2(l-e' 2 ) 1/2 (1 + e' 2 )} cos(2E*- 2g - 2h) 

+ e' {(2 - e' 2 ) - 2(1 - e' 2 ) 1/2 } cos(3E'- 2g- 2h) 

- 15 e' {(2 + e' 2 ) + 2(1 - e' 2 ) 1/2 } cos(E'+ 2g + 2h) 

+ 3 {(2 + e' 2 ) +2(1 -e' 2 ) 1/2 (1 +e' 2 )} cos(2E'+ 2g + 2h) 

- e' {(2-e' 2 ) + 2(1 - e' 2 ) 1/2 } cos(3E'+ 2g + 2h)]} . 

If AL = L - L' 

AG = G - G' 

AH = H -H', 
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the short periodic perturbations in the Keplerian elements a, e, I are 
obtained by 


I n * 

a = a + 2a 


al 

L' 


t 

e = e 




AG\ 
G # / 


Irl'+COtl' 


/ag_ah\ 

\ G f K # / 


and in the expressions for AL , AG , AH, the variables / , g, h may be 
replaced by the variables /' , g' , h' with an error of the 4th order. 

It is important to remember the fact that the primed variables are 
only affected by long periodic perturbations (with an error of the 4th 
order). 


12. The Long Period Terms. Elimination of h' 


At this stage the Hamiltonian F' = Fq + F' + depends on the 
variables g # , h' , L' , G' , H' , and 1/ is a constant with respect to time. 
The next step in von Zeipel's method consists of the elimination of h' 
and g' . The elimination of h' is performed by means of a generating 
function 


s' =l "i 9 + G"g' + H #/ h # +s; + s; + ••• 

of a canonical transformation from ( / ', g' , h # , L' , G # , H' ) to ( l " , g" , 
h" , L w , Q" , H" ). The new Hamiltonian 

F" = Fq + F'j + F" + ••• 


should be independent of h w . 
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In order to simplify the formulas the following symbols are intro 
duced: 


i = — = /l - e' 2 
L' 


G' 


Then 


F' = 


V- 


2L 


'2 


n?L' 4 f 

n;H' + -1(5 - 3 t7* 2 ) [(-1 + 30' 2 ) 

16 m 2 e l 


+ 3(1 -8' 2 ) cos 2h'] + 15(1 -V 2 ) 


1(1 +6') 2 cos2(g' + h') 


+ (1 - 6 >' 2 ) cos 2g' + 1(1- 9') 2 cos 2(g 


’ ' h > ]} 


^ [- b! Jj (1 - 3e ,J ) + 6 b 2 J„ (1 - 6' 2 ) cos 2h‘] 

4 l'6 \G / 


The elimination of h' depends on the solution of the system 


Fq(L") = Fq(L'') = * 


2L" 2 

f;(H') = FJ(H') =n;H‘ 
bs: BF,' 


f; = f; + • 


Bh' BH' 


and since Sj does not depend on l' , L 1 


= L\ 


As before a particular solution is given by 


f ; = 


F' 

r 2s 


= part of F' independent of h' 


bs: 

n^ = -F; p = -(F'-F' s ). 


Then it is found that 


F" - 
r 2 ~ 


n c L 


1 4 


16 ^ ? 


{(5 - 3t/' 2 )(-1 + 36 " 2 ) 


where 


+ 15(1 - 7?" 2 ) (1 - 0" 2 ) 


V 


cos (1-39-2), 


// 


CF 

L' 


6" 


ir 

G" 


Furthermore, 

n ' s; = -^£ ($ (I (S - 37, ’ 2) (1 - e ' 2) sin 2h ' 


+ ^ (1 - V 2 ) [(1 + 9 “) 2 sin 2 (g' + h') - (1 - 9 ") 2 sin 2 (g' - h')] j 


“ 4 n 2 b 2 J. 


22 


,"3 


(1 - d" 2 ) sin 2 h' 
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The long periodic perturbations depending on the motion of the node (and 
eventually on that of the pericenter) are given by 


c3s; 

T ' _ T " , 1 - T " 

L - L + B/' " L 


bs; 

G* = G" + r~ 7 = 

7g 


f 15 

/ -?\ 

1— 

CN 

V 

1 

i-H 

( “ 32e 

W n '/ 

II 

77 L 


ll i ll 


bs: 


H - H + 7—7 = 


Bh' 


f, 3 

l n * 'l 

1 r 

[ " 16 e 

Uc n 7 

77" 6>" L 


+ f C 1 - V" 2 ) {(1 + 6") 2 cos 2 (g" + h") 


t (1 - e ") 1 cos 2(g- - h')}] -|j£ (■£) b J J 22 i2^pcos 2h* 


i - e " 2 


v 


bs; 

/'=/•- 1 


bl" 


= / + 


1 I n ■ 


16£ Vn*n'' 


1 3(10 - 3 rj " 2 ) (1 - 6>" 2 ) sin 2h'‘ 


+ 15 ^1 j [C 1 + ^”) 2 sin 2 (g" + h") - (1 - 0") 2 sin 2(g" - h")] 


9n^ K2T 1 - g" 2 . 

>■ * b J 00 ho * i sin 2h 



bs; 

BG" 


= g + 


16e 



3 ( 50" 2 - 3t/' 2 ) sin 2h" 


-fa+e") (v" 2 +0") Sin2(g" + h") 


+ -y (1 - O W J - O sin 2(g" -h')j- 


3 n 


* 2 


4n* 


b 2 J 22 


3 - 5 0" 2 
G" rf 3 


sin 2h" 


h' 


= h" - 


bs; 

BH" 


= h" 




-5) sin 2h" 


+ — (1 - rj ' 2 ) [(1 + 6 ") sin 2 (g" +h") + (1 - 9 ") sin 2(g" - h") 


-|^b 2 J 

2 n* 


22 


g " 2 

H" T )" 3 


sin 2h". 
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It is important to note that these terms are factored by 1st order 
factors. Then, a second order theory produces first order long period 
terms. It is necessary to go to 3rd order to obtain second order long 
period terms. This is done next. 


13. The Second Order Long Period Terms and Elimination 
of h' and the Time 


Care must be exercised in the evaluation of the third order part of 
the Hamiltonian. One of the dangers involved is including in the third 
order part of the Hamiltonian terms which are really second order. 

This difficulty could arise, for example, in the evaluation of the solar 
perturbations. This would result in the introduction of more terms into 
the second order Hamiltonian, which is permissible. 

An additional complication is the uncertainty of the values of the 
zonal harmonics for the potential field of the moon. There is no harm 
done in the theory if J 2 and J 22 are really smaller than the values 
assigned to them. However, if J 4 is introduced and treated as though 
it were of order J* when it is really greater than J 2 , then factors of 
the form J 4 /J 2 might invalidate the theory. 

Keeping this in mind, one can list the small parameters of 3rd order 

(moon's potential field) 

(earth's oblateness) 


(sun's perturbation) 


(correction due to the inclination of 
the moon's orbit to its equator) 

(solar radiation) 
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(eccentricity of the moon's orbit) 



(earth's perturbations - 3rd 
Legendre polynomial) 


— 2 (correction due to physical libration) 

n 


Since short periodic perturbations of ttle third order will be neglected, 
it is understood that all terms of third order which depend on l will be 
ignored. Suppose that F 3 is the third order part of the Hamiltonian and 
,F 3 ' the part free of short period terms. The Hamiltonian will be time 
dependent through the longitudes of the earth and of the sun. This fact 
introduces one more degree of freedom in the problem and it is neces- 
sary to introduce the time as a canonical variable. This is done through 
use of the following equations (r = t + const.): 


dr _ 3? 

dt 3T 


so that 


and 


dT = 33; = 33; _ ^3 j 

dt + 3 t 3t 3t 

3 r = F - T 
T = - F 3 + const . 


This additional complication does not affect the development of the 
theory up to this point, if one sets 


3S' 


3t 


= 0 . 


55 


The next step is then to eliminate h # to 3rd order to obtain 2nd 
order long period perturbations depending on the argument h' . 

The orders of magnitude are as follows: 


5 ’ = Fi - T =-^7 - T 
0 0 2L 


=F; =n*H' 


V - F' 
J 2 “ r 2 


= f; • 


where T has been incorporated into for reasons that will become 
obvious . 

The von Zeipel differential equation of the 3rd order is then 


as; ay; as ;ay; as; a?; as; a?; 
3 + ah' ar + a g ' dG" + ah' aH" + ar at 



as; a?; 
acr V” 


» 


or 


as;ar; ar;as; * as; as; 
3 + ah' aH" + aG" a g ' + ah' ~ a r 


aF! as; 
- f" + — - — - 
3 a g ' aG" 
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Since we had supposed that BS 2 /Br = 0 in the previous order evaluation, 
♦ a particular solution of the 3rd order equation is found by taking 


/Bs; BF'\ 
F 3 = F 3s + [W W'l 


UG" Bg ' / 


bf" bs; 

Bg' BG" 


where the subscript n s n designates non-dependence upon h" and r. 
Further, 


~ * 
41 <c 




dH ’7 P h' 


/bf^ bs;\ /be; bs;\ 
\B G " Bg ' / ph < + \Bg ' B G "/ ph ' 

and 

BS' /Bs; BF,'\ /bf' bs;\ 

bT =F 3p- + BH 7 ' ,/ pT \BG" bT 7 ]^ 


M K\ 

\ Bg ' BG'/, 


where the subscript ph 7 indicates the inclusion of all trigonometric 
terms with arguments of the form jh # + kg', j and k integers (j^O), 
and the subscript pr indicates the inclusion of all trigonometric terms 
with arguments of the form j r + kh ' + mg' , where j, k, and m are 
integers ( j ^0) . 

Then will be obtained as a periodic function independent of the 
time, and as a periodic function dependent on time. 
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Since the functions F 2p , Fj , S [ , and F" are known, the only parts 
that have to be computed are those corresponding to the 3rd order terms 
of the Hamiltonian (F 3 ) . 


14. Third Order Terms Generated by Coupling 
of 2nd Order Terms 

It was found that 


F' = _L n'L' >1 (5 - 3v' 2 ) [(- 1 + 3 6' 2 ) + 3(1 - 6' 2 ) cos 2 h'] 

1 16e \n / 1 


+ 15(1 - 77 ' 2 ) 1 (1 + 6' ) 2 cos 2(g' +h') + (l-0' J ) cos2g' 


4 / t\ 3 

+ cos2(g- -h')]} [- b 2 J, (1 - 30' 1 ) 


+ 6b 2 J 22 (1 - 6' 2) cos 2h ' ] 


{(5- 37,-2) (_1 + 3 0"2) 

16e \n / 


+ 15(1 - v " 2 ) (1 - e" 2 ) cos 2g"}-± b 2 J, 


(1 - 36" 2 ) 



and 


_ n *s; = - _ 

® 1 16 



^(5- 377 " 2 ) (1 -d" 2 ) sin 2h ' 


+ — (l-v" 2 ) 1(1 + 8") 2 sin 2( g ' + h')-(l-£") 2 sin 2(g ' - h ' )] 


- n ' 2 b2 J 22 “l ( 1 - " 2 ) sin 2 h' , 

and in the coupling terms the variables in F 2 can be double -primed 
with an error that is of the 4th order at least. 

The partials needed are given below. 

w- - it. n ' (£)’ Vb- {2(S - ^ e ' 2 (1 - cos 2h ’> 

+ 5(1 - 77 " 2 )6" [(1 + 6") cos 2(g" + h") - 26" cos 2 g" 

- (1-6") cos 2(g" -h")]> + | b 2 n ' 2 [J 2 - 2J 22 cos 2h"] 

A H 77 


BF' 


2 3n ' f n z \ 2 1 


T^" = TT (-') “T (1(V" 2 -56" 2 ) + 2(S6" 2 - 377 " 2 ) cos 2h" 
<3G 16e \n / rj 


- 5(1 + 6") ( 77" 2 + 0) cos 2(g" +h") + 10 (6" 2 - t?" 2 ) cos 2 g" 


+ 5(1 -6") (6" - 77 " 2 ) cos 2(g" - h" )} (J 2 (1 - 56" 2 ) 

4 G "^«3 

- 2J 22 (3 - 5 6" 2 ) cos 2h".} 
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dS i i 3L ' / n l 


3h* - ^ V-32I {2(5 " 3V } (1 ‘ * } C ° S 2h " 


+ 5(1 - V' 2 ) [(1 + 6") 2 cos 2 (g" + h") + (1 - 6 ") 2 cos 2 (g" - h")]} 


3 n ' 2 u2 T 1 - 6 " 2 

- — b 2 J 22 cos 2h 

2 n l V'3 


is ( n l \ 1-7 j" 2 


3g' 32e \n*ny 77 " 


G" [(1 + 0 ") 2 cos 2(g" + h") 


- (1 - 0") 2 cos 2 (g" -h")] 


°F 2 15 / n \ 1 

^---Te n,V \t > ) “ 77 " 2 ) (1 - 0 ” 2 ) sin 2 g" 


**;_ 3 hi) 1 

3G" " 32e \n*n'/ 77 " 


(2 (50" 2 - 3 77" 2 ) sin 2 h" 


- 5(1+ 6")(t" 2 + ( 9 ")sin 2 (g" + h") + 5(1 - 0") ( V " 2 - 6") sin 2 (g" -h") } 
o n * ^ b 2 T 

_3 ~22_ (_3 + 50" 2 )sin 2h" . 

4 n * 77 " 3 G " 

From these equations it can be shown that 


S' (coupling) = + 


/ \2 / \2 

(h rS\ \L { _ 32 0 " (1 -0" 2 )t7" 2 (15-17t 7" 2 ) sin2h" 

4096 e 2 \n 7 In* / 77 " 


- 40" (1 - 0" 2 ) (5 - 3 t 7" 2 ) 2 sin 4h" 


80(1 + 0") 2 (2-30")t7" 2 (1 - 77 " 2 )sin( 2 g"+ 2h") 
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+80(1 -6>") 2 (2 + 30")t ? " 2 (1 -t/' 2 ) sin(2g" -2h") 

+10(1 + d") 2 (1 - 6")(l-r]" 2 ) [5 <l-d")-6r," 2 ] sin(2g"+4h") 
+ 10(1 -0") 2 (1 +i9")(l -T)" 2 ) [5(1 +i9")-67/' 2 ] sin(2g" - 4h") 
+25(1 +i9") 3 (1-t?" 2 )[(1-i9")-2t?" 2 ] sin (4g" +4h") 
+25(1-^") 3 (1-7 7 " 2 )[(1+^")-2 t 7 " 2 ] sin(4g" - 4h" ) } 

9 /n \ 2 

+ _|_ n '7/'- 4 Hn J 2 b 2 {40"(l-0" 2 )(5-3r ? " 2 )sin2h" 

+5(1+<9") 2 (1 +2<9"- 5<9" 2 )(1- 7/' 2 )sin(2g" + 2h") 


+5(1-9") 2 (1 - 29" - 59" 2 )(l-7 ? " 2 )sin(2g" - 2h")} 

g / n \ 2 

r— n H Hr J 2 2 b 2 {-89"(1-6>" 2 )(5-3 t/' 2 ) [2 sin 2h" - sin 4h"] - 

51 2 e \n*J 

-20(1 + i9") 2 (1 - 8")(3 - 59") (1 -t/' 2 ) sin (2g" +2h") 
-20(1 -6") 2 a +9")(3 +59")(1 -t]" 2 ) sin (2g" - 2h") 


+5(1 +9") 2 (1 - 9")(1 + 59")(1 - 7]" 2 )sin(2g" + 4h") 


+5(1 -<9") 2 (1 +9") (1 -50") (1 -7]" 2 ) sin (2 g" - 4h") } 


4 ( 4 ] -•* 


8 n* 


16 K 



« 6 jj ,; 
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a' 2 /l-e" 2 cosI"{ [(2+33 e" 2 )-(2-l 7e" 2 )cos 2 I-"] 


F" (coupling) 




rr 


+ 1 5e" 2 sin 2 I"cos 2g" } 


+ 


_9 

32 


/ n <t\ / n «\ n ' 2 sin 2 1" cos I 
' n 'VW e(l - e" 2 ) 2 


J 22 b 2 {2(2 + 3 e" 2 )+15 e" 2 cos 2 g" } 


9 /n'\ n' 2 sin 2 I" cos I" 
+ 4 [ n *J a' 2 (l -e" 2 ) 7/2 


and, o£ course, 


S' (coupling) = 0. 


It is important to note that (coupling) is made up of second order 
terms though it is obtained by means of multiplying a second order quan- 
tity by a first order quantity. 

The next few sections are devoted to the computation of F 3 ' . From 
that computation various* additional parts of , S 3 , and S' are then 
derived. 


15. The Radiation Pressure 


For this computation the shadow effect will be neglected. Its in- 
clusion would require, of course, the introduction of a shadow function 
due both to the moon and the earth. 

Consider a to be the absolute value of the acceleration of the or- 
biter arising from the solar radiation. Then, the disturbing function 
for the radiation pressure will be given by 


R a = - ar cos S" 3 
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or, neglecting the inclination of the lunar equator to the sun's orbit, 


R -~crr[cos(f + g)cos(ft- X ) 


- sin (f + g) sin(ft - X ) cos I] 


The elimination of short periodic terms gives 

=-| era' e' jeos 2 J-cos(g' +h' + iV 0 - ^ @ ) + sin 2 h cos (g' - h' + \ 0 - K m ) |» 


where use was made of the formulas 


_L f 

2 ” 1 


2tt 


r cos f d/ = - — a e 

2 


2tt 

r sin f d / = 0 , 


and where the node ft was written in terms of the canonical vari- 
able h. Thereisno secular contribution from this effect, that is, 

F^ s (radiation) = R^ s = 0. On the other hand, since there are no terms 
strictly independent of time, the contribution to is zero, or 



S' (radiation) = 0. 


Therefore , 


BS^ (radiation) 

Br 



(radiation) = , 


or 


S' (radiation) 


3 a n r j ti 

crael o J- • /// i// \ \ * 

+ ^-* co S 2 sm (g +h 

2 n t- n tt 2 

- sin 2 ^ sin (g" - h" + - ^ e ) j ’ 
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where 


= n ® r + const • 
x. e = n * r + const . , 

and n* and n* are, respectively, the mean motion in longitude of the sun 
around the moon and of the moon around the earth. Of course, n*-n* ~ n 
is a good approximation. 


16. The Second Legendre Polynomial for the Sun’s 
Gravitational Perturbations 


To restate what has been said previously, the most important term 
in the disturbing function due to the sun’s gravity is 

t, (£) ,p , 

where, for the computation of S" 3 , the sun can be considered to move 
along the moon's equator. 

Then 

cosSj 3 = cos (f + g) cos (fi - \ 0 ) - s in (f + g)sin(fi-\ 0 )cosl. 


On the other hand, neglecting the mass of the moon, 


n 


2 

© 


k 2 (m Q + m e ) = k 2 m 0 



and since the ratio m 0 /m Q can be neglected, 


n ® a 3 = ^3 
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0 * 


Furthermore, if the eccentricity of the earth orbit is neglected, 


F 3 (sun) = n 2 r 2 P 2 (cos S" 3 ) = 

= ^n 2 r 2 (3 cos 2 S" 3 - 1) = 

= ^n 2 r 2 {cos 2 (f + g) cos 2 (Q -iV 0 ) - 2 sin (f+g) cos (f + g). 

• sin(fi-X. @ )cos (Q- X. 0 ) cos I + sin 2 (f+g) sin 2 (Q-\ 0 ) cos 2 i} 



That this contribution is of the third order is visible from the fact that 



where 

^0 ^ 2 S _ 
n £ ~ 365 

Now 

(sun) = ~ n 2 r 2 {(cos 2 f cos 2 g + sin 2 f sin 2 g) cos 2 (ft - A. 0 ) 

- 2 (- sin 2 f sin gcosg + cos 2 f sin g cos g) sin (H- \ Q ) cos (ft- A. 0 ) cos I 
+ (sin 2 f cos 2 g + cos 2 f sin 2 g) sin 2 (ft - X. 0 ) cos 2 I } 

- I n ® 1-2 + f n ® r2 sin f cos f * Q * 

where Q is independent of f. 
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In order to obtain the part of F 3 (sun) independent of /, F 3 (sun) is aver 
aged over /. The following well known integrals are used: 




27 T 


2 


r‘ cos 


f d / = a 2 ( — + 2 e 2 


n 


2 77 


r 2 sin f cos" f d/ = 0 


(n = 0 , 1 , 2 , ... ) 


r.[ 


27 T 


r 2 sin 2 f d/ = — a 2 ( 1 - e 2 ) 

2 


rj 

" r 


2 77 


r 2 d / = a 2 


( i+ f 


They are obtained using the eccentric anomaly as the variable of inte- 
gration. It follows that 


F 3 ' (sun) = — n 2 a " 2 | - 2 (2 + 3e" 2 ) (1 - 3 cos 2 I") + 30 e " 2 sin 2 I" cos 2 g" 
j z 


+ 6(2 + 3e" 2 ) sin 2 I" cos 2 (h" + \ 9 -k Q ') 


+ 15e" 2 (1 + cos I ") 2 cos 2 (g" + h" + _ ^©) 


r "\2 


+ 15e" 2 (1- cos I ") 2 cos 2 (g" - h" - V 0 + A. @ ) j 


The secular contribution is 


F; (sun) = — n 2 a " 2 {-(2 + 3e" 2 ) (1-3 cos 2 I") + 15e" 2 (1 - cos 2 I") cos 2g 



There is no contribution to Sj.and the contribution to S' 3 is 
S 3 < sun > = ^ 3T—T ( 2 ( 2 + 3e2 > sin ' l ' si n 2 (h* ■+ \ 0 - \ e ) 

n (E n © 

+ 5e" 2 (1 + cos I") 2 sin 2 (h" + g" + - \ @ ) 

+ 5e" 2 (1 - cos I") 2 sin 2 (h" - g" + - \ @ )} 


17. The Third Legendre Polynomial for the Earth' s 
Gravitational Perturbations 


Referring to page (19) 


Now, 


F 3 (earth) = * P 3 (cos S/ 0 ) . 

n 2 a| = k 2 (m e + m e ) = k 2 (l e, 

so that considering the moon's orbit to be circular, 
F 3 (earth) =^^P 3 (cos S/ 0 ) = 


C .3 


e a. 


r P 3 (cos S/ 0 ) 


1 

2 


n: 


e a„ 


(5 


cos 


s ; 0 - 


cos 


Si'o), 
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where it is permissible to write 


cos Sj 0 = cos (f+ g) cos (fi- A. e ) - sin (f + g) sin cos I - 

= -^cos f [(1+cos I) cos (g+h) + (1-cos I) cos (g-h)] 
- i sin f £(l + cos I) sin (g+h) + (1-cos I) sin (g-h)] . 


If 



+ cos I) cos (g+h) + (1-cos I) cos 


(i-h)] 


then 


so that 


B 0 = i [(1+ cos I) sin (g+h) + (1- cos I) sin (g-h)] 


cos Sjq = A 0 cos f - B 0 sin f, 


2 

F 3 (earth) = \ — — r 3 {5 [A 3 cos 3 f- 3A 2 B 0 cos 2 f sin f 

1 e a <t '• 

+ 3A 0 B 2 cos f sin 2 f-B 3 sin 3 f] - 3 [A 0 cos f-B 0 sin f]| 


The short period terms are eliminated with the aid of the integrals 


f.J, 


2 rr 


- 3 cos 3 f dl = - 


/ 1 5 5 3 \ 

V 8 6 + 2 6 / 


- — a 3 e (3 + 4 e 2 ) 


i r 


r sin 


in f cos" f d/ = 0 


(n = 0, 1, ... ) 
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ro| H- 


It follows that 



cos f sin 2 f d/ = - — a 3 e (1 - e 2 ') 

8 v ' 


cos f d/ = - -| a 3 e (4 + 3 e 2 ) . 
8 


F 3 (earth) = - — 


n l a' 3 
e a„ 


{s[(3 t 4e' 2 )A;*t3AiB; 2 <l-e' 2 )] 


-3A; (4+3e’ 2 )}, 


where 

A ® 3 = 32 ^ “ 6 ' 2 ^ ( 1-0, ) cos (g'-3h') + (1-0') 3 cos (3g'-3h') 

+ 3'(1- 0' 2 ) (1-0') cos (3g'-h') + 3(1-0') (3+ 20' + 30' 2 ) cos (g'-h') 
+ 3 (1+0') (3- 20' + 30' 2 ) cos (g' + h') + 3(1-0' 2 ) (1+0') cos (3g' + h') 
+ 3 (1-0' 2 ) (1 + 0') cos (g' + 3h') + (1 + 0') 3 cos (3g' + 3h')} 


(1+0') cos (g' + h') +I(l-£») cos (g'_ h ') 


K B; 2 = j- 2 {(1 + 0') (3- 20' + 3d' 2 ) cos (g'+h') 


+ (1-0') (3+ 20' + 30' 2 ) cos (g'-h') - 3(1-0') (1-0' 2 ) cos (3g'-h') 


+ (1-0') (1-0' 2 ) cos (g'-3h') - (1-0') 3 cos (3g'-3h') 
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- (1+0 ') 3 cos (3g' + 3h') - 3(1 +e')(l-6' 2 ) cos (3g' + h') 
+ (l+e')(l-d' 2 ) cos (g' + 3h')}. 


The contributions to F 3 and to S 3 are both zero. The contribution to 
S' 2 is then given by 


Si (earth) = 




(earth) dh' , 


or 


S' (earth) - 1536 ( 


5 _/S\2i 


„»3 

a e 


\ n c/ 


a c e 


{-9(l + cos I")(l + 10cos I"- 15 cos 2 I")(4+3e" 2 )sin(g"+h") 


+ 9 (l-cos I" )( 1-10 cos I" -15 cos 2 I")( 4 + 3 e" 2 )sin(g"-h") 
+ 15 sin 2 I" (1+cos I")(4+3e" 2 ) sin (g" +3h") 

-15 sin 2 I" (1-cos I")(4+3e" 2 ) sin(g"-3h") 

+ 315 sin 2 I* (1 + cos I") e " 2 sin ( 3 g" +h") 


-315 sin 2 I" (1-cos I") e " 2 sin (3g"-h") 


+ 35 (1+cos I ") 3 e " 2 sin (3g" +3h") 
-35 (1-cos I ") 3 e " 2 sin (3g"j3h")} . 


18. The Eccentricity of the Moon’s Orbit 


The correction for the eccentricity of the moon's orbit is given by 


(e.) 




r 2 P 2 (cos S 10 ), 
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or, keeping only the first power of e £ , 



r 2 cos / 

© 


P 2 (cos S 10 ) 


3 e c n c 

4 € 


r 2 (3 cos 2 S 10 - 1) cos / 9 


where the mean anomaly l 9 is = n £ t + const. Then, 

2 *2 

F 3 ' (e c ) = ^ "* ** e e {- 2 (2+3e' 2 )(l-3 cos 2 I') + 30e' 2 sin 2 I' cos 2g' 
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+ 6 (2+3e' 2 ) sin 2 I' cos 2h' + 15e' 2 (1+cos I # ) 2 cos (2g' + 2h') 


+ 15 e' 2 (1- cos I') 2 cos (2g'-2h')} cos,/ ( 


There are no contributions to S' or F". The contribution to S' 3 is 


3 n <c a e <c 


^ - 64 4 


{- 2 (2+3 e" 2 ) (1-3 cos 2 I") 


+ 30e ,/2 sin 2 I"cos2g" + 6 (2+3e" 2 ) sin 2 I"cos 2h" 
+ 15 e"? (1+cos I") 2 cos (2g"+2h") 

+ 15 e" 2 (1- cos I"-) 2 cos (2g"-2h")| sin / 0 . 


It is worthwhile to note that although this result is a part of a third 
order generating function its actual order is (n^/n) e c which is consid- 
ered to be second order. A small divisor n c is introduced through the 
integration. 
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19. The Inclination of the Moon's Orbit to its Equator 


The expression for the earth's perturbation is 


n e 1-2 


P 2 (cos S' 0 ), 


where cos S ! 0 = K(i ) • By a Taylor expansion, if sin i e is small, 


, ^K(i_) , 


= cos S' 10 + 


Boos S! 


10 


cos F B 
<c c 


sm 


*e-° 


Thus the correction to be introduced is 


3n£ r 2 B cos S! 0 

i. w. F sin cosS '» 


d c 


‘c - 0 


From the expression for the cos S' Q it is easily found that the quantity 
in the bracket is 


sin I sin v sin (f + g). 


Thus 


n c r a 


P 3 C i(gl = 3 cos s ; 0 sin (f + g) sin v e sin I sin 


Now 


cos S' Q sin (f + g) = sin ( f + g) cos ( f + g) cos h 


- sin 2 (f + g) sin h cos I 
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Since 


i_ [ 2n 

277 Jo 


r 2 sin (f + g) cos (f + g) d l ~ — a 2 e 2 sin 2g 


2v I 


277- 


r 2 sin 2 (f + g) d / 


-ii r 

2 l 2 " J„ 


r 2 d / 


2 - Jn 


r 2 cos 2(f + g)d/ 


= -j { a2 i 1 + f e2 ) " f a2 ® 2 cos 2g }’ 


the part of F 3 (i £ ) which is free from short periodic terms is 


. 3n 2 a' 2 

^3 = — v ® sin I' sin i e { - 2(2 + 3e' 2 ) cos I' sin h' 


8e 


+ 5 e' 2 (1 + cos I') sin(2g' +h') + 5e' 2 (l - cos I') sin(2g' - h') } 


There is no contribution to F" and S' • The contribution to S' is 

o Z 3 

3n 2 a" 2 

S 3 (i^) = - ~ — cos v 0 sin I' sin i £ { - 2(2 + 3e" 2 ) cos I" sin h" 


8e(n c + Naj £ ) 


+ 5e ,/2 (l + cos I") sin(2g w + h") + 5e" 2 (1 - cos I") sin(2g" - h") } . 


where N co^ is given by a) = t + const. /Again, this is a second 

order contribution since n 2 / (n. + NoO ~ n^. 

<£ v Z T z z 
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20. The Non-Sphericity of the Potential Field of the Earth 


Because of the fact that the earth's equator is not the reference 
plane, the form of the disturbing function due to the zonal harmonic co- 
efficient j 2 of the earth will be derived from basic relations. 

In Figure 5, E is the earth, M the moon and S the orbiter. 


S 



Figure 5 


The disturbing function for the motion of S is given by 



P 2 (sin cjp) 


where R 0 is the equatorial radius of the earth and cp the latitude of the 
orbiter with respect to the earth's equator. The part k 2 m 0 /r / has al- 
ready been taken into account and the rest is supposed to be a 3rd order 
quantity. Therefore 


f 3 (.) = . 



Rj j 2 Pj (sin <p>. 


If the terms (a/a c )j 2 and e^ j 2 are neglected, it follows that 

n £ R J j 2 

F 3 (®> " _ l P 2 ( Sin ^ ’ 

assuming a € /r' 3 = 1. 
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The angle 9 must now be expressed in terms of the orbital elements 
» of the lunar equator. 


In Figure 6 the geometry of the problem is given. 



Figure 6 

F rom basic relations of spherical trigonometry it is found that 
sin cp = (sin I cos I 0 - cos i p cos I sin I $ ) sin (f + g) 

- sin I e sin <// cos (f + g) . 

Then, if 

A = sin I cos I 0 - cos \p cos I sin I ffi 
B = _ sin I 0 sin i p, 
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we have 


f 3 (®) = - 




{[j <A 2 + B2 > - \] * 1 < B2 * A2 > 


cos 2 (f + g) 


+ | AB sin 2 (f + g) 

In order to eliminate the short periodic terms, the following integrals 
are needed: 



where 


2tt 


1_ 

2tt 


i 


2n 


cos 2 fid/ -2 / 3 2 


2 TT 

\ sin 

Jo 


2fd / = 0, 



/ 3 = - (l - /l - e 2 ). 
e 

(The use of/3 here is not to be confused with its use earlier as the latitude.) 
They are easily obtained by considering the expansion 

00 

d/ = df + 2 ^ k/3 k ^ + /l - e 2 j cos kf (-l) k df, 
k=l 


which is the differential of the equation of the center. It follows that 


+ |^ 2 (i * /^) 

• [ (B 2 _ A 2 ) cos 2g' + 2 AB sin 2g' ]|- » 


where 


y3' = /3(e'). 
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I 


Furthermore, since \f> = h + \ 0 - Cl', where fi' is the longitude of the 
•descending node of the moon's equator on the earth's equator, we have 

B 2 + A 2 = ^sin 2 I' cos 2 I 0 + | cos 2 I' sin 2 I 0 + I sin 2 lj\ 


1 

2 


sin 2 I 0 sin 2 I' cos 2 (h' + 


fi') 



sin 2I 0 sin 21' cos (h' 


fi') 


B 2 _ A 2 = sin 2 I 0 - sin 2 I' cos 2 I 0 -_L cos 2 I' sin 2 I 0 

\2 2 


-i- sin 2 I 0 (1 + cos 2 I') cos 2 (h' + - fi') 


+ -j sin 2 I' sin 2 I 0 cos (h' + k @ - Cl'), 

and 


AB = -| sin 2 I e cos I'sin2(h' + \ 0 _ Cl') -1 sin2 I 0 sinl' sin (h' + \ 0 - fi') . 
Thus, the contribution to is zero, and the contribution to is 



+ 3/3" 2 (l +2 /l - e" 2 ) sin 2 I" cos 2g"} (1-3 cos 2 I 0 ) . 
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The contribution to S ' 3 will be given by the integration of 
BS' (®) 3 n l h R l sin 1 


Bt 


16 


| _ 2 sin I" [sin I" sin I 0 cos 2 (h' + ) 

+ 4 cos I" cos cos (h' + ^ ) ] 


+ /3" 2 (l + 2 A - e" 2 ) [-(1 + cos I") 2 sin I 0 cos 2(h' + g' + A-© - ft* ) 

- (1 - cos I"') 2 sin I 0 cos 2(h' - g' + k 9 - ^ ) 

+ 4sin I" (1 + cos I") cos I 0 cos(h +2g' + k® _ Q' ) 

- 4sin I" (1- cosl")cosl 0 cos(h' - 2g' +^ 0 - ft') ] }• 


Since 


we have 


Bt 


<\» . 0 .') = n* - N n 


c ’ 


= - _L n|i ^ sinI » 
32 «< n I - N n«' 


2 sin I" [sin I" sin I 0 sin 2(h"+ \ 0 - ft' ) 

+ 8cos I" cos I 0 sin (h"+ k “ ft ) ^ 


/ 3 " 2 (l + 2 / 1 - e" 2 ) [-(1 + cos I") 2 sin I sin 2 (h" + g" + - ft') 


-(l-cosI") 2 sin I sin 2 (h"- g" + k - 


+ 8sinI"(l + cos I") cos I 0 sin(h" + 2 g" + k# - ft' ) 
-8sinI"(l-cos I" - ) cos I 0 sin(h"- 2g" + - ft')] j • 
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I 


2 1 • Higher Order Zonal Harmonics for th e 
Moon’s Potential Field 


The terms corresponding to J 3 , J 4 , J 5 are obtained immediately 
from Brouwer’s theory for artificial satellites (Reference 1). Their 
contribution is entirely to the "secular" part of F' . So 

O n' 2 J 3 b 3 

F 3 (J 3 , 4 . 5 ) = - 4 e"sin I" (1 - 5 cos 2 I") sin g" 

8 a'(l _ e " 2 ) 5/2 

3 n ' 2 J 4 b4 

— [(3 -30 cos 2 I" + 35 cos 4 I")(2 + 3 e" 2 ) 

128 a ' 2 (l -e" 2 ) 772 

- 10(1- 8 cos 2 1" + 7cos 4 I")e" 2 cos 2g" ] 

5 n ' 2 J s b5 

- 7 TF 7 e" sin I" [6(1 - 14 cos 2 I" + 21 cos 4 I") • 

256 a' 3 (l _ e" 2 ) 972 

• (4+ 3e" 2 ) sin g" - 7(1 -10 cos 2 I" + 9 cos 4 I") e" 2 sin 3 g" ] . 


The order s of magnitude of these terms will depend, for reasonable val- 
ues of e", on the values of J 3 , J 4 , and J 5 , which are not very well 
known. 


j 

22. Physical Libration, and the Precession of the 
Lunar Equator 


Physical libration causes a periodic oscillation in the position of 
the lunar surface, thereby creating a small angular displacement be- 
tween the principal lunar meridian and the earth-moon line of centers. 
The largest contribution to this displacement is given by 

X = a sin / t ( _ ) + x 0 (see Reference 9, p. 316) 
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where 


a = - 59" = _ 2.86 xl0- 4 rad. 

X 0 -constant dependent upon the initial time 
- mean anomaly of the sun. 


Let f 77 £ and xyz be lunicentric, equatorial coordinate systems 
with the £ -ax i s directed toward the earth and the x-axis passing through 
the principal meridian (see Figure 7). 


Then, 



Figure 7 


" _ BF _ BF 


x = £ cos y + 7] sin x 
y = - £ sin y + 77 cos X- 


Thus , 

x = £ cos X + V- si * 1 X + 2 an 0 cos -£ 0 (77 cos y - £ sin y) 

-an* sin -^(77 cos y - £ sin y) - a 2 n| cos 2 -{. 0 (t 7 sin y + £ cos y), 
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where n @ = mean motion in the mean anomaly of the sun. 


V cos x - £ sin x = an w * cos l + y 


•• . or 

<f cos x + 7] sin x = , 

Bx 


+ 2an 0 cos (an @ x cos -i 0 + y) - an 2 sin (77 cos x - f sin 

- a 2 n 2 x cos 2 


~O(10“ 12 ) 


~O(10- 9 ) 

n 2 


^ ~ O(10- 6 ). 


Therefore, if one neglects terms higher than third order, one has 


•• BF _ - 

x 1 = — + 2an y cos -i 
Bx 


A similar computation yields 


y = 2 an 0 x cos ■i Q - 


Thus , 


.. BF .. _ BF 
Bx By 
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where 


F = F + 2 an 0 (xy - xy) cos -l Q 

= F + 2 a n 0 H cos . 

Hence, the addition to the Hamiltonian of the problem is 


F 3 (lib. ) = 2 an Q H cos . 

Then, 

F 3s (lib.) = 0 
F 3pl/ (lib.)*0 
F 3 P r^ lib ^ : 2an 0 H cos -i Q , 

and therefore 
BS' (lib. ) 

— = 2an 0 H cos S' 3 (lib. ) - 2a H sin -£ 0 . 

Now, the precession in the lunar equator is created by the small 
angle of inclination ( ~1 °32*. 1 ) between the lunar equatorial plane and 
the ecliptic. This results in the regression of the equatorial node in 
the ecliptic. However, the motion of the node is already implicit in 
the n* _ hence, no additional corrections need be made. 


23. Complete "Secular" Third Order Hamiltonian 

The complete secular third order Hamiltonian is given by 

F 3 = F" (coupling) + F 3 (sun) + F 3 (®) + F" (J 3>4i5 )- 

Formulas for the various parts of F" can be found in preceding 
sections . 
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24. The Determining Function for Long Period Terms 


Additional long periodic perturbations depending on the motion of 
the node, the motion of the earth and the motion of the sun will be given 
through the determining functions S' and S' . The determining functions 
S' and S^ are 


S' = S' (coupling) + S' (earth) 

S' = S 3 (radiation) + S' (sun) + S' (e^) 


+ S 3 (i e ) + S 3 (©) + S' (libration) . 


In the next section, the partial derivatives needed in order to find 
additional perturbations in the canonical elements are given. 


25. Long Period Perturbations of Second Order 


The perturbations of long-period are obtained from 


l 


II 


BS 

3L" 


= / 


bs; ^ ^ 

BL" + BL" + BL" ’ 



= L" 


g 


II 



!5 + Z!a bs^ ^ 
BG" BG" + BG" ’ 


h" = 


BS 

BH" 


h' 




+ * 


bh" bh" 




= G" + 


= H" + 


bs; bs' bs; 

Bg' + Bg' + Bg' 


bs; 

BS' 

bs; 

Bh' 

+ Bh' 

+ Bh' 


The terms corresponding to Sj have already been obtained. Next 
the partial derivatives of S' and S 3 with respect to a" , e", I", g' , 
and h' are computed. For convenience, however, the primes have been 
dropped in this section. 
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The "mean” mean motion n' (called n) depends on a' (called a) 
through the relation 


n 


2 



and 


3S _ 2 /a 3S | 1-e 2 BS 

BL* ' ^ Ba e Be 


BS _ 1 /l - e 2 BS t cot I BS 

BG e V Be / ^,a ( 1 — e^) ai 


BS _ _ 1 BS 

sin i / /+a (1 - e 2 ) ^1 


The partial derivatives of the various parts of and S^ are as follows: 


BS^ (coupling) 63 

= + 


8192 e 2 


1^.] (A\ Q£_ { 32 sin 2 I cos I (l-e 2 )(2 

V n / \n*J /HP 


- 17 e 2 ) sin 2h 


- 4 sin 2 Icos I (2 + 3 e 2 ) 2 sin 4h 


+ 80(l+cos I) 2 (2-3 cos I) e 2 (1-e 2 ) sin (2g + 2h) 

+ 80 (1-cos I) 2 (2 + 3 cos I) e 2 (1-e 2 ) sin (2g-2h) 

+ 10 sin 2 I (1+cos I) e 2 [6e 2 - (l + 5cos I) ] sin(2g + 4h) 
+ 10 sin 2 I (1-cos I) e 2 [6e 2 -(l-S cos I) ] sin (2 g-4h) 
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+ 25 (1+ cos i) 3 e 2 [2e 2 ~(l + cos I)] sin (4g+4h) 

+ 25 (1- cos I) 3 e 2 [2e 2 -(l~ cos I)] sin (4g-4h) } 


j b 2 2 

512 € V n c/ 2 a(l-e 2 ) 2 


{4s in 2 I cos I (2 + 3 e 2 ) sin 2h 


+ 5(l+cosI) 2 (l + 2cosI-5cos 2 I)e 2 sin (2g + 2h) 


+ 5(l-cos I) 2 (1-2 cos I - 5 cos 2 I)e 2 sin (2g-2h) } 


27 / \ 2 n sin 2 I 

1024 e \n*/ J « a(1 _ e2) 2 


{ -8 cos I (2 + 3e 2 )(2 sin 2h-sin 4h) 


-20(l+cos I)(3-5 cos I) e 2 sin (2g + 2h) 


-20(l-cos I)(3 + 5 cos I) e 2 sin (2g-2h) 


+ 5(l+cos I)(l + 5cos I)e 2 sin (2g+4h) 


+ 5(1- cos I)(l-5 cos I) e 2 sin (2g-4h) } 


117 

16 



n sin 2 I cos I 
a 3 (l-e 2 ) 7/2 


sin 2h 


117 / n \ 2 ^2 b 4 n s ^ n2 ^ cos I 
32 {<) 22 a 3 (l-e 2 ) 7/2 


sin 4h 
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BSj (coupling) 

31 

• {- 16(1-3 cos 2 I)(l-e 2 )(2-17e 2 ) sin2h 
+ 2(1-3 cos 2 I)(2 + 3e 2 ) 2 sin 4h 

- 40 (1 + cos I)(l-9 cos I) e 2 (1-e 2 ) sin (2g+2h) 
+ 40 (1-cos I)(l +9 cos I) e 2 (1-e 2 ) sin (2g-2h) 
+ 10(1 + cos I) e 2 t (3- 5 cos I)(l + 2cos I) 

-3e 2 (l-3cos I)] sin(2g+4h) 

- 10 (1-cos I)e 2 [(3+5cos I)(l-2cos I) 

-3e 2 (l+3cos I) ] sin(2g - 4h) 



+ 25 (1+ cos I) 2 e 2 [2(1+ cos I) -3e 2 ] sin (4g + 4h) 

- 25 (1- cos I) 2 e 2 [2(l-cos I) -3e 2 ] sin (4g-4h) } 


9 J 9 b 2 nsin ^ {-(1-3 cos 2 I)(2 + 3e 2 ) sin 2h 


64 e W/ 2 ( i . e2) 2 

- 5(l+cos I )( 1-cos I-5cos 2 I)e 2 sin (2g+2h) 

+ 5 (1- cos I)(l+ cos I- 5 cos 2 I) e 2 sin (2g-2h) } 
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00 I o 


9 

256 e 



BS' (coupling) 


n \2 . j 

— 1 J, 2 b 2 n Sin - { 4(1-3 cos 2 I)(2 + 3e2)(2 sin 2h-sin 4h) 
n c/ (1-e 2 ) 2 

- 20(l+cos I)(l + 7 cos I~10cos 2 I) e 2 sin (2g + 2h) 

+ 20(1- cos I)(l-7 cos I-10cos 2 I) e 2 sin (2g-2h) 
-5(l+cos I)(l + 2 cos I)(3-5 cos I) e 2 sin (2g+4h) 


+ 5(l-cos I)(l-2cos I)(3+5 cos I) e 2 sin (2g-4h)} 


4 n 55 T (1-3 cos 2 I) ou 

J 2 J 2 2 b 77Z SAn A11 

* 2/i 2\7/2 

a z (1-0 


t2 t 4 n sin I (1-3 cos 2 I) . A , 
J 2 b 4 1 L sin 4h 

22 2 /i 2 n7/2 

a z (1-0 


J 


+ 9 (M 2 /M 2 n 

4096 e 2 ' \ n t) ( 1 


2 

e 


(l-e 2 ) 3/2 

•{-96 sin 2 I cos I (12-17e 2 )(l-e 2 ) sin 2h 


-4 sin 2 I cos I (2+3e 2 )(14-9e 2 ) sin 4h 
+ 80(1+ cos I) 2 (2-3 cos I)(2-3e 2 ) (1 - e 2 ) sin (2g + 2h) 

+ 80(1- cos I) 2 (2+3 cos I)(2-3e 2 )(l-e 2 ) sin (2g-2h) 

+ 10 sin 2 1(1+ cos I) [6e 2 (4-3e 2 )-(2-e 2 )(l + 5cos I)] sin(2g+4h) 
+ 10 sin 2 1(1- cos I) [6e 2 (4-3e 2 )-(2-e 2 )(l-5 cos I)] sin (2g-4h) 
+ 25(1 + cos I) 3 [2e 2 (4-3e 2 )-(2-e 2 )(l+ cos I)] sin(4g+4h) 

+ 25(l-cos I) 3 [2e 2 (4-3e 2 )-(2-e 2 )(l- cos I)] sin(4g-4h)} 
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J 2 b 2 


ne 


1 28 e V n*J J2 (W)' 


{ 4 sin 2 I cos I (7 +3e 2 ) sin 2h 


+ 5(l + cos I) 2 (l + 2cos I-5cos 2 I)(l+e 2 ) s in (2g+2h) 


+ 5( 1-cos I) 2 (l-2cos I-5cos 2 I)(l+e 2 ) sin (2g-2h)} 


— 2_ j b 2 ne sin2 ~ {-8 cos I (7 + 3e 2 )(2 sin2h-sin4h) 

256e \n l) 22 (l- e 2 ) 3 


-20(1+ cos I)(3-5cos I)(l + e 2 ) sin(2g + 2h) 


-20(1- cos I)(3+ 5 cos I)(l+e 2 ) sin (2g-2h) 
+ 5(1+ cos I)(l + 5 cos I)(l + e 2 ) sin (2g+4h) 


+ 5(1- cos I )( 1-5 cos I)(l + e 2 ) sin (2g-4h)} 




ne sin 2 I cos I 


a 2 (l-e 2 ) 


2 . 9/2 


sin 2h 


. 63 / n \ 2 J2 b 4 ne sin 2 I cosl sin 4h 
16 Ki “ a 2 (l-e 2 ) 9/2 


SSj (coupling) 

3g 


+ 16 — V /— V n a2 e2 { 8(l + cos I) 2 (2-3cos I)(l-e 2 ) 

1024 e 2 V n ) [ n *J v/fT? 


• cos (2g+2h) 


+ 8(1- cos I) 2 (2 + 3cos I)(l-e 2 ) cos (2g-2h) 

+ sin 2 I (1+cos I) [6e 2 -(l + 5cos I)] cos (2g + 4h) 


+ sin 2 I (1-cos I) [6e 2 -(l-5cosI)l cos (2g-4h) 


+ 5(1+ cos I) 3 [2e 2 -(l + cosI)] cos (4g+4h) 


3S' (coupling) 
3h 


+ 5(l-cos I) 3 [2e 2 -(l-cos I)] cos (4g-4h) } 


45 /M 2 
128 e In* ' J2 


J 2 b 2 — 2-2 — { ( 1 + cos I) 2 (1 + 2 cos I-5cos 2 I)cos(2g + 2h) 
«/ (1-e 2 ) 2 


+ (l-cos I) 2 (l-2cos I-5cos 2 I)cos (2g-2h)} 


45 / n «V T , 2 ne2 sin 2 I , ... T ... _ T v 

J, 9 b 2 — {-4(1+ cos I)(3-5 cos I) cos (2g+2h) 


256 e \n* J J22 


(1-e 2 ) 2 


-4(1- cos I)(3 + 5 cos I) cos (2g-2h) 


+ ( 1 + cos I)(l + 5 cos I) cos (2g+4h) 


-l- (1- cos I)(l-S cos I)cos (2g-4h) } 



{16 sin 2 Icos I (l-e 2 )(2-17e 2 )cos 2h 


-4 sin 2 I cos T (2 + 3e 2 ) 2 cos 4h 


+40 (1+cos I) 2 (2-3cos I)e 2 (l-e 2 )cos(2g+2h) 


-40( 1-cos I) 2 (2 + 3 cos I) e 2 (1-e 2 ) cos (2g-2h) 


+ 10 sin 2 I (1 + cos I) e 2 [6e 2 -(l+5 cos I) ] cos (2g + 4h) 


- 1 0 sin 2 I (1-cos I) e 2 [6e 2 -(l-5cosI)] cos (2g-4h) 


+ 25(l+cos I) 3 e 2 [2e 2 - (1 + cos I) ] cos (4g + 4h) 


-25(l-cosI) 3 e 2 [2e 2 -(l-cos I) ] cos (4g-4h)} 
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— i — (— ) J b 2 {4 sin 2 I cos I (2 + 3e 2 ) cos 2h 

128 e\n* / 2 (1 _ e 2 } 2 


<Ci 


5(1+ cos I) 2 (l + 2cos I-5cos 2 I) e 2 cos (2g+2h) 


-5(1- cos I) 2 (1-2 cos 1-5 cos 2 I) e 2 cos (2g-2h)} 






n sin 2 I 


1 28 e \n*/ (1-e 2 ) 2 


{-8 cos I (2 + 3e 2 )(cos 2h-cos 4h) 


-10(l+cos I)(3-5 cos I) e 2 cos (2g+2h) 


+ 10 ( 1-cos I)(3+5cos I) e 2 cos (2g-2h) 


+ 5(l+cosI)(l + 5cos I) e 2 cos (2g+4h) 


-5 (1- cos I)(l-5 cos I) e 2 cos (2g-4h) } 


9 /n \ 2 t , 4 nsin 2 Icos I 
+ - (-*)J 2 J22 b — cos 2h 


4 \n 




a 2 (l-e 2 ) 7/2 


9 / n \ 2 t2 , 4 n sin 2 I cos I 


4 In 


J22 b 4 


2/i 2 \ 7/2 

a 1 (l-e^) 


cos 4h 


BS; (earth) 5 / n c\ n c a ‘ e 


Ba 


512 e\n $J a. 


{-9(l + cos I)(l + 10 cos 1-15 cos 2 I)(4+3e 2 ) sin(g+h) 


+ 9(l-cos I)(l-10cos 1-15 cos 2 I)(4 + 3e 2 )sin (g-h) 


+ 15 sin 2 I (1 + cos I)(4 + 3e 2 ) sin (g+3h) 



-15 sin 2 I (1-cos I)(4+3e 2 ) sin (g- 3 h) 


+ 315 sin 2 I (1 + cos I) e 2 sin (3g+h) 
- 315 sin 2 I (1-cos I) e 2 sin (3g-h) 


+ 35(1+ cos I ) 3 e 2 sin (3g+3h) 


- 35 (1-cos I ) 3 e 2 sin (3g-3h) } 


(earth) 5 /n £ \n £ a 3 e sin I 

zi = 5127 


(3(ll-10cos I-45cos 2 I)(4+3e 2 )sin (g+h) 


+ 3(11 + 10 cos I-45cos 2 I)(4 + 3e 2 ) sin (g-h) 


-5(l+cosI)(l-3cos I)(4+3e 2 ) sin (g+ 3h) 


-5(l-cosI)(l+3cos I)(4+3e 2 ) sin (g_3h) 


- 105(l + cos I)(l-3 cos I) e 2 sin (3 g+h) 


-105(l-cos I)(l + 3 cos I) e 2 sin (3 g-h) 


-35 (1 + cos I ) 2 e 2 sin (3g+3h) 


-35(1- cos I ) 2 e 2 sin (3g- 3h) } 


3S' (earth) 5 /n £ \ n £ 


51 2 e \n*J a £ 


{-3(l+cos I)(l + 10 cos 1-15 cos 2 I)(4+9e 2 ) sin (g+h) 


+ 3(l-cos I)(l-10 cosI-15cos 2 I)(4+9e 2 ) sin (g- h) 


+ 5 sin 2 I (1+ cos I)( 4 + 9 e 2 ) sin (g+3h) 
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- 5 sin 2 I (1- cos I)(4+9e 2 ) sin (g-3h) 


+ 315 sin 2 I (1+cos I) e 2 sin (3g+h) 


- 31 5 sin 2 I ( 1 -cos I) e 2 sin (3 g-h) 


+ 35 (1+ cos I ) 3 e 2 sin (3g + 3h) 


- 35 (1- cos I ) 3 e 2 sin (3g-3h) } 


Bs 2 (earth y _ — 5 _ /n^\ n c a f {_ 3 (1 + C os I)(l + 10cos I -15 cos 2 I)(4 + 3 e 2 ) cos (g+h) 
3g 51 2 e \n*y a c 


3 (1- cos I)(l-10 cos 1-15 cos 2 I)(4+3e 2 ) cos (g-h) 


+ 5 sin 2 I (1+ cos I)(4+3e 2 ) cos (g+3h) 


- 5 sin 2 I (1- cos I)(4+3e 2 ) cos (g- 3 h) 


+ 315 sin 2 I (1+cos I) e 2 cos (3g+h) 


-315 sin 2 I(l-cos I) e 2 cos (3 g-h) 


1 35(1+ cos I ) 3 e 2 cos (3g + 3h) 


-35(l-cosI) 3 e 2 cos (3g-3h) } 


BS' (earth) c /n.\n a 3 e 

— - = — [— !•) — {-3(1+ cos I)(l + 10 cos 1-15 cos 2 I)(4 + 3e 2 ) cos (g+h) 

Bh 512 e\n*/ a c 

-3(1- cos I )( 1-10 cos 1-15 cos 2 I)(4+3e 2 ) cos (g-h) 


+ 15 sin 2 I (1+ cos I)(4+3e 2 ) cos (g+3h) 
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+ 15 sin 2 I (1-cos I)(4 + 3e 2 ) cos (g-3h) 
+ 105 sin 2 I (1+ cos I) e 2 cos (3g+h) 

+ 105 sin 2 I (1- cos I) e 2 cos (3g-h) 

+ 35(1+ cos I) 3 e 2 cos (3g + 3h) 

+ 35(l-cos I) 3 e 2 cos (3g-3h) } 


“S 3 (radiation) 3 ae f j 


Ba 


= I n *_n* 1 C0S2 f Skl (H+g+ + Sin ^ t Sin ( h -g + A -®-^®)} 


C ® 


BS' (radiation) 3 crae sin I 


BI 


4 "*-n: 


{- sin (h+g+ \ e -X. @ ) + sin (h-g+ \ 9 - \ e ) } 


BSj (radiat ion) 3 


Be 


= "2 n : a n , {cos 2 ! sin (h+g+\ e -\ 0 ) + sin 2 | sin (h-g+ \ e -\ 0 )J- 


BS' (radiation) 3 




= | {cos 2 I cos (h+g+\ $ -\ e )- sin 2 !cos (h-g+\ e -\ 0 ) j 


^S' (radiat ion) 3 

3h 


= 2 { C ° s2 \ C ° S (h+g+ X ®~ k + Sin2 7 COS (h ~ g+ ^©)} 


c ® 
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(sun) 

3a 


BSjCsun) 

31 


3S' (sun) 

3e 


3S^ (sun) 

3g 


3S^ (sun) 

3h 


— — (2(2 + 3e 2 ) sin 2 I sin 2 (h + - \ 0 ) 

32 < - n ® 

+ 5e 2 (l +cos I) 2 sin2(h + g + \ 9 - >^ Q ) 
+ 5e 2 (l - cos I) 2 sin 2(h - g + /V e - 


o n? a 2 sin I 

— { 2(2 + 3e 2 ) cos I sin 2 (h + k - 

32 

- 5e 2 (1 + cos I) sin 2(h + g + X 0 


+ 5e 2 (1 - cos I) sin 2 (h - g + \ 


3 ^ a e , 

— {6 sin 2 I sin 2 (h + - k ) 

32 <-< 

+ 5(1 + cos I) 2 sin 2(h + g + \ e - A- 
+ 5(1 - cos I) 2 sin 2(h - g + - A 


15 

32 


n: 



{(1 + cos I) 2 cos 2(h + g +\ 0 - A.©) 


_ 3 _ 

32 


nia 2 


* 

n <L~ n e 


- (1 - cos I) 2 cos 2(h - g + k e - X 


(2(2 + 3e 2 ) sin 2 I cos 2 (h + k® - A 


+ 5e 2 (1 + cos I) 2 cos 2(h + g + k^ 


+ 5e 2 (1 - cos I) 2 cos 2(h- g + \ 


BSi (O ^ 

-I- n s e s a {-2(2 + 3e 2 )(l - 3 cos 2 I) + 30 e 2 sin 2 I cos 2g 


Ba 32 e 


+6(2 + 3e 2 ) sin 2 Icos2h + 15e 2 (l + cosI) 2 cos (2g+ 2h) 


+ 15 e 2 (1- cos I) 2 cos (2 g - 2h)} sin l m 


BS; (e 4 ) 9 


Bl 32 e 


n $ a 2 sin I {- 2(2 + 3e 2 ) cos I + 10 e 2 cos I cos 2g 


+ 2(2 + 3e 2 ) cos Icos2h-5e 2 (l+cosI)cos(2g + 2h) 


+ 5e 2 (l-cosI)cos(2g- 2h) } sin L 


9 

Be 32e 


n e c a 2 e {-2(1 - 3cos 2 I) + 10 sin 2 Icos 2g 


+ 6sin 2 Icos2h + 5(1 + cos I) 2 cos(2g + 2h) 


+ 5(1 -cos I) 2 cos(2g - 2h)} sin L 


BS 3 < e <t ) 45 


Bg 


32 e 


n £ e a 2 e 2 {2sin 2 I sin 2g + (1 + cos I) 2 sin(2g + 2h) 


+ (1 -cos I) 2 sin(2g - 2h) } sin ^ 


dS 3 9 


Bh 


32 e 


n (E e (E a2 {2(2 +3e 2 )sin 2 I sin 2h + 5 e 2 (l +cos I) 2 sin(2g +2h) 


- 5e 2 (1 - cos I) 2 sin (2g - 2h) } sin L 
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96 


3s; (®) 

3a 


3S' (®) 

31 


3S' (®) 
3e 


3s; (®) 


o 


3 “c 

16e 


Jo R £ sin L 


n^-N 


(1<C 


{-2[sinIcosI sin I 0 sin 2 (h + \ 0 - fl' ) 


-4(1-2 cos 2 I) cos I 0 sin(h + A_ 0 - Cl* )] 

+ yS 2 (1 +2/ 1 -e 2 ) [sin 1(1 + cos I) sin I 0 sin 2 (h + g + \ 0 - Q' ) 

- sin 1(1 __ cos I)sin I 0 sin 2(h - g + \ 0 - fi') 

-4(l + cos I)(l-2cosI)cos I 0 sin(h+2g + \ 0 - fl' ) 
-4(1- cos I) (1 + 2 cos I)cos I 0 sin(h - 2 g + A. 0 - ft')] } 


3 n ®j 2 R ®^ 2 sin 1^(2 +/l-e 2 ) 

{-Cl +cos I) 2 sin I sin 2(h+g+\ - D') 

66 

- (1-cos I) 2 sin I ffi sin 2 (h-g+A. 0 - O') 

+8sin 1(1 +cos I) cos I 0 sin(h+2g + \ 0 - F?) 
-8sinl(l-cosl)cosl 0 sin(h-2g + ^ 0 - fi')) 


o n l J 2 KP 2 sinI ® C 1 + 2/1-e 2 ) 

{-(l+cosI) 2 sinI cos2(h + g+\ - ^') 

6e e ( n £ - Nq^ ) 

+(1 -cos I) 2 sin I 0 cos 2(h- g + \ 0 - fl') 
+8sinI(l-|-cosI)cos I 0 cos (h + 2g + \ 0 - .0') 
+8sinI(l-cosI)cos I 0 cos (h-2g + - fi') } 
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BS; (e) 

Bh~ 


_ _3_ n I J 2 R 2 0 sinI e { _ 2 sinl [sin I sin I 0 cos 2(h + \ 0 - ft') 

166 n «“ Nqj 

+ 4cos I cos I 0 cos (h + - O' )] 

+ / 3 2 (1 + 2 /l-e 2 ) [-(1 +cos I ) 2 sin I 0 cos 2 (h + g + A. e - G') 


-( 1 - cos I") 2 sin I 0 cos 2 (h - g + X 0 - fl ) 

+ 4sinl(l + cos I) cos I 0 cos(h + 2g + \ 0 - ) 

_ 4 sin I (1 - cos I) cos I 0 cos (h - 2g + ^ ) 


as; cub.) =0 

Ba 


BSl (lib. ) . T . p 

3 = _ 2 aGsinI sini 0 

31 


BS; (lib.) 

_1_ = 0 


BS; (lib.) 

3g 


= 0 


BSl (lib.) 

_ 0 

Bh 


26. Secular Perturbations and Perturbations Depending Strictly on g" 
At this stage of the problem the Hamiltonian is 

r = % + t; + y; + y 3 - 
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, where 


= — - - T = const. 

0 2L' 2 

= n*H" = const. 

?2 = *£(-, g", - , L', G" , H" ) 

= F 3 . g", V , G" , H") 

and L' , T and H" are constants in time. T is effectively a 3rd order 
quantity but in the manipulation of von Zeipel’s method its contributions 
are of the zero order, for it is a variable itself, independent of all the 


It is easy to see that the variable part of the Hamiltonian 3" , which 
will generate the equations of motion, is factored by small parameters. 
Therefore, a method of successive approximations such as von Zeipel's 
cannot be applied. 


Therefore, the system of differential equations produced by 3" 
must be integrated directly or by using some kind of approximation. 

The 3rd order part, 3" , will be neglected for the moment. The remain- 
ing Hamiltonian is 

F" =3" +FJ + F" = const. ' 

where, since L' , H" and T are constants, F" = const. The equations 
of motion are 


3F" 

G = = G (g", G") 

°g 


g" =~^r=k"U", G") 


h = ■ n ® ■ w = h (g ''’ G,,) 
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1 = n * - Tl j=l (g", G") 


The first two equations represent a system of one degree of freedom 
since F" is a first integral. Then it can be integrated and therefore 
h" and /" are obtained by quadrature. 

The approach taken in this section is similar to the one used by 
Kozai in Reference 6 . The equations for if and g" are found to be 


*)"= -5 K x (l - V" 2 ) (1 - 9 " 2 ) sin 2g" 


g" = — i [( 77" 2 - 5tf" 2 ) - 5 ( 77 "" - 0"") cos 2g"] - 


,/2 


L t i/4 

V 


(1 - 5 0 


and the integral 6F" is written 


KjL' [(5 - 3 V " 2 ) (-1 + 3d" 2 ) + 15 (1 - v" 2 ) (1 - 0 " 2 ) cos 2g" 


(1 - 3 d" 2 ) = C = const. = 6F" 


v 3 , / n <t \ 2 3 , 1 /a' \ 3 

K. = — — n I — = — n a. . a. =— — 


1 ~8e 4e 1 ’ 1 ~ 2 a. 


T7- _ 3 1 2 , 2 T 3 / t 2 

K - — n tr J = - n L a cl 

2 4 12 4 2 * 2 ,2 

a 
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In order to integrate the system, the following remarks are made: 

(1) 77" is assumed to be a periodic function of time of the 

form 


17" = 77" + A77" (a v a 2 ; t ) , 

where Vq is a constant and A77" strictly periodic in t and such 
that A77" ( 0 , 0 ; t) = 0 for any _t , so that 

A77" = a l A77" (a r a 2 ; t) + a 2 A77" (a,, a 2 ; t) • 

(2) g is assumed to be the sum of a linear function of time plus 
periodic terms, i.a.. of the form 

g" = g" + n g (a v a 2 )t + Ag" (ctj, a 2 ; t) , 

where g<, is a constant, n g is a constant which can eventually 
be zero (libration cases) and Ag" is strictly periodic in t . They 
are such that n ( 0 , 0 ) = 0 and Ag" ( 0 , 0 ; t) = 0 for any t . 

Thus 


n g = a l n i ( a i» a 2^ + a 2 n 2 ( a i> a 2^ 

^g — a j Agj ( a ji a 2’ ^ + a 2 ^§2 ^ a i ’ a 2’ ^ * 


( 3 ) If assumption ( 1 ) is correct, we have 

V = a 1 At7 x + a 2 At7 2 , 

and, comparing with the equation for 77", it is seen that A 7 7" 
must be zero, or say At 7"(0, a 2 ; t) = 0 for any a 2 and t. Then 
77" will depend on <x 2 through a fourth order term of the form 
a i a 2 ' ’ * ° ne concludes that if a 2 is neglected in the evaluation 
of 77" , the error will possibly be a periodic function factored by 
the product 0^ a 2 , and it is tolerable. 

The same is not true for g" since neglecting terms ctj a 2 * • 
may not be tolerable due to the presence of a linear term in 
time. 
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27. The Integration of rj 


Rewriting the integral F" we have 

(5 - 3t/ /2 X-1 + 3 6" 2 ) + 15(1 - rj" 2 )(l - 6 " 2 ) cos 2 g" = 

= 1 [C + 2K t ]"- 3 (1-3 e" 2 ) ] = W . 

KjL' 2 

If 7 j“ is the value of rf for g” = 0 , then neglecting the oblateness 
terms we have 


where 


W = 10 - 12 rj' 2 + 6 v 2 = — - — = const. , 

1 KjL' 


v - — = cons t. 

v 


We are going to evaluate cos 2g" from this integral in order to find 
the value of dx/dt. The error in doing this is of the fourth order by 
remark (3) in the previous section. Therefore, 


cos 2g" 


2(5 - 6xj + 3v 2 ) + (5 - 3x) (1 - 3 v 2 x -1 ) 
15(1 - xXl - v 2 x _1 ) 


where x = rj" 2 and x x = rj'^ 2 . Note that the definition of r | x ' 2 does not nec- 
essarily correspond to a physical situation. 


The substitution of cos 2g" into the equations for rj' gives 



Also 


AjG" 



BF 2 

— dt = L" (/x - v^ 0 ). 


In the equation for dx/dt the minus sign is used when sin 2 g" is 
positive and the plus sign is used when sin 2 g" is negative. The choice 
of sign will be treated in the next section. 

The values of x 2 and x 3 are defined by the relations 


■^2 "t ^3 — + 5 V 2 — 2 Xj ) 

and, for eccentricity e" < 1 , they are always real. 

The solution for x will be generally expressed by elliptic functions, 
except in a few cases where x x , x 2 , x 3 assume special values. 

The possible cases for the relative values of x x , x 2 , x 3 are next 
determined based on the fact that one looks for a real solution of the 
problem. Any complex solution, from the physical point of view, must 
be discarded. 

It is noted that, since 


x = rf 2 = 1 - e" 2 

= i!ll = (i _ e " 2 ) cos 2 I", 
L' 2 

x and v satisfy the inequalities 

1 - x ^ v 2 - 0 . 
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Furthermore, if the solution x has to be real, then from the differential 
equation for x there is the additional condition 

(x - Xj) (x - x 2 ) (x - x 3 ) £ 0 • 


( a ) x t = 1 
Then 

If 

then 

and since x £ 1, the 



x i = x 2 = X 3 = 

nly possibility is 

x = 1 = const. 


or 


e" = const. = 0 


cos I" = ± /VS = const. 


If 


then 


2 3 - 

v z >— , X, > 1, 

5 J 

1 = Xj = x 2 < x 3 , 


and again the only possibility is 

x = 1 = const. 


or 


e" = const. = 0 


cos I 


> / 3/5 • 



then 


1 - Xj - x 2 > x 3 

and 


Thus , 


or 


x 3 Sx ^ 1. 

— y 2 < x £ 1 

3 

0 < e" < 1 


0 ^ cos I" ^ 1. 


(b) Xj = v 2 
Then 



and 


0 ^ x 3 = x x = v 2 ^ 1 . 


The only possibility is 


or 


x = v 2 = const, 
cos I" = 1 and/or e" = 1. 


(c) v 2 < Xj f 1 
Thus , 

v 2 < x 2 < , x 2 ^l 

and 

^ 2 <x 3 <|, x 3 / 

This is the most general case. 
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28. General Case ( v 2 < x x / 1) 


In this case, Xj, x 2 , x 3 will be ordered in ascending order and will 
be called x', x' 2 , x' 3 so that x 3 > x' 2 > x' • 

The solution for x is expressed in terms of Jacobi's elliptic 
functions as 


x - X 1 + ( X 2 “ X j) 


snucnvdnvt snvcn udnu 


1 - k 2 sn 2 u sn 2 v 


where 
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and cnv, dn v , sn v are the Jacobi's elliptic functions with modulus k 
and argument v ; x Q is the value of x at the intial time which is taken 
equal to zero . 

Since at this stage g" is unknown, the choice of the sign can be 
made only after the equation for g" is integrated. That is, the correct 
sign is determined by a comparison between the value of the integral 
F" using cos 2 g" obtained from the integration and its value at t = t Q . 
The solution for x given above is one in which the oblateness portion of 
the energy integral has been neglected. By an iteration process, the 
correction due to oblateness can be taken into account. 


29. The Integration of g” 


Once x is a known function of time, g" can be integrated using x 
as the independent variable. In fact, consider the transformation 


dg; 

dt 


dg" dx 
dx dt 


dg" 

= + 4 /6 Kj /Q(x) > 

dx 


where 


Q (x) = (x - Xj) (x - x 2 ) (x - x 3 ) . 


On the other hand, 


• // 


g 


K i K 0 

[ (x - 5^ 2 x* 1 ) - 5 (x - v 2 x' 1 ) cos 2g" ] (1 - 5v 2 x" 1 ) 

/x L' X 2 
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and 


cos 2g" = - (• — ^ + 2 

15 (1-x) (1-v 2 x- 1 ) |K 1 L K,L' 


2 v-l 


2 1 - 3v x 

v 3/2 


+ (5 - 3x) (1 - 3 v 2 x ~ 1 ) V » 


so that 


g" = -K 1 /x+5K 1 i' 2 x -3/2 


C x - v 2 x -1 


3L' /x (l-x)(l-v 2 x _1 ) 


^ ^2 (x-v 2 x~ 1 )(1-3v 2 x *) 

3L 7 X 2 (1-x) (1- V 2 x -1 ) 


(x-v 2 x ')(5-3x)(l-3v 2 x *) 
3 /x (1 - x) (1 - V 2 x _1 ) 


K 2 1 - 5v 2 x _1 

V 7 2 


After reduction, this becomes 


• ;/ 


g 


-£./ i ^ — — 1 ah / ( i.3v»)-i 

3L' \(l-x) /x (x-v 2 )/x /x J ^ L |_ 



7 + 6 v 2 1 | 9 v 2 \ 
2 x 2 + 2 x 3 J 
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or, 


+ — J 2 (1 - 3 ^ 2 ) 2v 2 (3v 2 -5) - 

3 l (l-x)/x (x-v 2 )/x 


- 2 (3v 2 + 1) 


— ) 

V^J 


finally , 


74 ^ K ^ = 3 [u t2(1 - 3 ^ K .] 


(1 -x) /xO(x) 


+ 2v 2 (3v2_ 5)Ki 


(x - v 2 ) /xQ (x) 




QOO 


- — |(i- 
3 L' t 


3v 2 ) L 


(1-x) /QOO y2 (x-v 2 )/Q00 


v 2 -3v 4 -2 1 7 + 6 


* 1 


* / QOO 2 x 2 /(foo 


9v‘ 


2 „3 


>=1. 

X 3 /Q (x) J 


109 


The integration gives 


+ 4/6Kj Aj g" + 2 (1 - 3 v 2 ) K x j 


dx 


, X() (1 -x) v/x Q (x) 


1 V Cv : 
“ 3 L L' 


+ 2 v 2 (3v 2 - 5)K, 


r 


dx 


(x-v 2 ) / x Q(x) 


-ire 

3 LL' 


+ 2 (1 + 3^ 2 ) Kj 


f X dx 
Jx n /x Q (x) 


*377 


r dx 4 «> r 

■L (l-x) / Q (x) 3L'v 2 J (x- 


dx 


v 2 )/QOO 


2 ^ 2 v 2 - 3 v* - 2 

+ 3 L 7 ^ 


•L x/Q(x) 3 L J x 


dx 

X 2 /Q (x) 


where 



dx 

, 

x 3 /Q(x) 



f (x)dx 



f (x)dx 



f (x)dx. 


The various elliptic integrals needed are given in Appendix A. 


30. The Variables /" and h" 


Once x is known, /" and h" are integrated in the same manner as 
was done for g” . The original equations are 
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{ (10- 3x) (- l + 3v 2 X -1 ) 


.N-^ = n ..n' 


31 / 


8e 



+ 15 (2 - x) (1 - v 2 x' 1 ) cos 2g" } - |Hl!b 2 J 2 (l-3^ 2 x- 1 ) x“ 3 / 2 

« L 


BF" 

BH" 


3 n' v 
8e x 



(5- 3x) 


- — n' 2 b 2 T — x _s/2 . 15 n* v ( n cY .. „ 

2 h V *71 ~ {-) (1 - x > cos26 . 


• „ / n .\2 


where cos 2g"is to be replaced by the formula in the previous section. 
Therefore 


l =n' - K . (10-3 x) (- 1 + 3 v 2 x -1 ) — — — — 

3 1 3L' 1-x 


2 ^2 (2-x) (1 - 3 v 2 x~ x ) _^i (2-x) (5 - 3 x) (1 - 3v 2 x -1 ) 

3 L* x^/2 (i _ x) 3 i^x 


^ 2 1 - 3 ^ 2 X -1 

L' 
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h" = - n * - K. 1/ 5 ~ — - 2K, — 

© 1 _ 2 T I 


,5/2 


Cl/ 1 2 K 2 v 1 - 3v 2 x- 1 

3L' x (i _ i/ 2 x _1 ) ^ L' x s/2 (i _ ^2 x -i) 


(5-3x)(l-3i/ 2 x _1 ) 

3 x (1 - i/ 2 X* 1 ) 

After reduction, these become 


n 


+ f K ‘ 


Cl 

3L'_ 


C 1 2 - 6v 2 

3L' 1 - x 3 1 - x 


2^2 | ~7-6i/ 2 1 1 - 3 V 2 2 1 1/ 2 1 

3 L* 2 x/x (1 - x)/x ^ x 2 /x” . 


h" = [-n; + 2 i/Kj] + 


l^v^i/ 2 - 5) + 


Ci/~ 

3L' 



2 r _2 _J_ _ _2 1 

3 L' _v 2 x /x v 2 (x - i/ 2 ) /x_ 

Finally, we have 

C 
3L' 


T4/6 - 

dx 


“ , 2 „ Cl 1 

n + — K- ■■ - 

3 3L'J /oTx) 
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K i 2 x 1 2 2 7 - 6 v 2 1 

* / r\ /- .1 \ __ . 

3 J (1 _ x ) A>(x) 3 L' L 2 x/xQ(x) 


1 - 3 v 2 21 v 2 1 

+ - ZT ~ r\ 

( 1 - x) v'xQ(x) 1 x 2 i/xQfx) . 


2 „2 


+ 4/6KJ ^1= (-n; + 2vKj)— 1=+ K, v(3v 2 - 5) + 

dx * /QOO L 3 3L - 


3Lj (x - v 2 ) v / Q(x) 


+ 2 V [2 1_ 


^ L' |_v 2 xv / xQ(x) v 2 (x - 


i—i 

! ) /xQ(x)J 


which, integrated, give 


+ 4/6 Kl A, = T4^K 1 (n' + |K 1 --^ 7 ) t- - 


K x (2 - 6 v 2 ) c 


(7-6 z^ 2 ) f — ^r + 7v 2 -J- f - 
J x /xQ(x) L J x . x 


/xQ(x) L J x x v xQ(x) 


- — — (1 - 3i^) f* 1— 

3 L' J (1 - x) i/xQ(x) 


r2K,v „ Cvl f 

: 4/6Kj Aj h" =+ 4/6K 1 (-n*+2K 1 v)t + |— — (3v 2 - 5) + — J J 


3vL> L J x (x - V 2 ) / xQ(x) Jx 0 X / xQ(x). 
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The various elliptic integrals needed are given in Appendix A. 

c 

3 1 . Special Cases 


The following special case remains to be treated: 



Since x x = x 2 > x 3 , the equation for x becomes 


dx 

dt 


= + 4 /6 K x 


(1 — x) Vx - X 3 


whose solution is 


x = x. + (1 - x.) 


1 + H exp [ + 4 \/~6 Kj ( — V 1 — x 3 t)] 
1 - E exp [ + 4 Kj ( - V 1 - x 3 t) ] 


where 


E = 


/x 0 " X 3 - 1/1 ' X 3 


/x 0 - x 3 + 1/1 - x 3 


and the choice of the sign — or + is done as in the general 
case after integration of the equation for g" . 
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The equations for g" , /" and h" are the same as for the general 
, case with the exception that in Q(x) , Xj = x 2 = 1 , so that the integra- 
tion is essentially different. In fact , /Q(x) = (1 _ x) v'x - x 3 , so that 


T4/6K 1 ^=+1 

1 dx 3 


-^2 + 2(1 - 3 v 2 ) K, 

Lt 


l 

(1 - x) 2 Vx(x - X 3 ) 


1 

3 


V 


+ 2v 2 (3v 2 -5) Kj 


1 

(x - v 2 ) (1 - x) Vx(x - X 3 ) 


1_ 

3 


+ 2(1 + 3 V 2 ) 

L 


i 

(1 - v'l Vxfx - x,~) 

' s o 


2 K 2 


1 2 
-+ ■ 


1 


+ - — J + (l -3v 2 ) 

^ L* L (1 - x) 2 v'x - x, v 2 (x - v 2 ) (1 _ x) v'x - X- 


- 3v 4 - 2 


1 


7 + 6u 2 


x(l — x) vx - x. 


2 „2 


(1 — x) Vx - X, 


9 v 2 


^ x 3 (1 - x) /x - x. 


T4/6K.^ = 

1 dx 


, 2 v C 

n ■(■ ■— K. — 

3 1 3L' 


J (1 - x) VX - X, 


■ c K, 

— + — (2 - 6i^ 2 ) 

3L' 3 


1 2 K 2 [V 

! v'x - x. 3 L' [ 


- 6 v 2 


(1 - x) 2 v'x - x 3 


'• x(l - x)Vx(x -X,) 


T (1 - 3v 2 ) 


21 v 2 


(1 - X ) 2 Vx(x - X.) ^ x 2 ( 1 


1 - 

— x) Vx(x - X 3 ) 
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3 L' 


v 1 2 (x - v 2 ) (1 _ x) /x(x - X.) 


The result of the integration is 


+ 4/6^ A, g" = + i ^ + 2(1-3v 2 )kJ X = 

11 3 L' JJ (1-x) 2 /x(x - x ) 

0 


1 Cl' 0 2 /o,i2 


3 L' 


+ 2v 2 (3v 2 - 5) K 


i £ - 


1 il + 2(1 + 3v 2 ) kJ f dX 

3 Ll' J Jx 0 (1 - x) /x(x - x 3 ) 


+~<1 > 


(1 - x) 2 vx - X 


(x-v 2 )(l - x)vx-x 3 °Lk 


. + |— (V 2 -3v 4 -2) j 

; 3 lv 2 J x(i-x) 


*2 . 7 + 6v 2 dx 

L' 3 *L x 2 (1 - x) /x-x. 


c 

L' 1 x 


(1 - x) VX - x 3 




The various integrals needed are given in Appendix B. 


This ends the particular cases to be considered. In resume, 

the precision of the results is as follows: 

(1) In cos 2g" to integrate x: error is proportional to a 2 

(Z) In the integrated value of x: error is proportional to 

a 1 a 2 > if no iteration is done to improve x. 

(3) In the integration of g", n , i: error is proportional 

to a]a 2 or a^a 2 . This means that the motion of the peri- 
center and node and the perturbations in the mean anomaly 
are, at maximum, in error by a sixth order quantity or 
about 10~ 12 . If the time elapsed is not too large, these 
errors are tolerable. 
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32. Higher Order Perturbations 


The solution of the problem up to this point stands as follows: 

(1) Secular terms: correct to 2nd order, included in the integra- 

tion of / ", g ", h". 

(2) Long period terms: correct to 2nd order in terms depending 

on h" and t (and eventually g" ). 

(3) Long period terms: correct to 1st order in terms depending 

strictly on g" and and ^ c . 

In order to evaluate long period perturbations correct to 2nd order 
in terms depending on g" , the 3rd order part of the Hamiltonian should 
be included. 

Consider then the equations 


V 


*3 

BL 


7= (e\ G") 


A 2 i" 


3F'; 


3F" 

A 2h" = -^=f 3 (g", G") 


The functions f. contain variables g" , G" factored by 3rd order param- 
eters (^10‘ 6 ). The variables g" and G"(x) depend on time through very 
complicated functions. The substitution of the independent variable in 
terms of x is out of the question, since t is a function of x through the 
inverse Jacobian function sn" 1 . The only possibility is then to consider 
G* (x) constant. If G" is a periodic function of time, the error in making 
this approximation is of the 5th order (10~ 10 ) in G" and could increase 
linearly with time in g" . Nevertheless, for a not too long time the error 
in g" may be tolerable since it is essentially 10~ 10 n't. For n' — 300°/ 
day, the error in longitude is less than 10“ 4 degrees for several hundred 
days. 
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I 


Therefore, if v" 0 , and g^ are the values of 77 ", 6" and g" at the 
epoch, it is possible to consider 

= g'o - |4 [ « - 5 > - 5 (V'o 2 - £'o 2 ) cos 2gJJ 

l^o 


K, 


L' < 


(l-5d''*)y t = g" + n" t 


G" = G'' 


If these values are substituted in f. ( i = 1,2,3), the 2nd order long per iod 
perturbations in l " , g", h" are obtained by quadratures, i.e. 


V" = / f, (g" (t), G") dt 


a 2 g " - / f 2 (g; (t), Gq> dt 


A 2 h" = f f 3 ( g ;' <t), G'')dt. 


For the variations in G" (3rd order), the equation is 

BF" 

A 2 G" =-^=§ (g", G") 
which can again be integrated considering 

A 2 G" .= /§(gj(t), G”) dt, 


where 

g; (t) = g'o + nj; t. 


Note that the quadratures above will give the secular perturbations to 
3rd order. 
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33 . Additional Long Period Perturbations of Second Order 
and Secular Perturbations of Third Order 


The third order part of the Hamiltonian, free from short periodic 
terms and from h" and r , is designated^ and appears in Section 22. 

In the partial derivatives that follow, the Keplerian elements are con- 
sidered double primed, but for convenience the primes have been dropped. 
The functions fj, f 2 , f 3 , and § are given next. 


BF 3 /a 3F 3 1-e 2 BF 3 


bf; i n 


- e" 

JJL3. 




cot I 




//it a (1 - e 2 ) 


BF" 

f 3 tg? (t). G''] - 3 


bf; 


s in I J~ix a (1 - e 2 ) 


g [g[ (t), G" ] = 


Bp" 


where 


= + 63 /V 
3a 2 56 l n* 


f n 


$ \ n 
n 


— a / 1 - e 2 cos I { (2 + 33 e 2 ) - (2 - 17 e 2 ) cos 2 I 
e 2 


+ 15 e 2 sin 2 I cos 2g } 


27 

64 



J 22 b 2 


ni sin * 1 c ? s 1 {2 (2 + 3 e 2 ) + 15 e 2 cos 2g> 
a (1-e 2 ) 2 £ 
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8 



n 2 sin 2 I cos I 
a 3 (1-e 2 ) 7/2 


+ — n 2 a {- (2 + 3 e 2 ) (1-3 cos 2 I) + 15 e 2 sin 2 I cos 2g } 

8 
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3 3 n 2 e sin I (1 - 5 co s 2 I ) ^ 

2 3 a 2 (l-e 2 ) s ' 2 


15 t ,4 n 2 


+ — — J* o' 

128 a 3 (1 - e 2 ) 7/2 


{ (2 + 3 e 2 ) (3-30 cos 2 I + 35 cos 4 I) 


- 10 e 2 sin 2 1(1-7 cos 2 I) cos 2g } 


15 

128 


n 2 e sin I 


+ J 5 b s -I ! — - ttlL-L— (6(4 + 3e 2 )(l-14cos 2 I + 21 cos 4 I) sin g 


a 4 (1 - e 2 ) 


2x9/2 


- 7 sin 2 1(1-9 cos 2 I) sin 3g } 


BF_" 

31 ‘ 


- ( — -\ (——] — a 2 / 1 - e 2 sin I { (2 + 33 e 2 ) - 3 (2 - 17 e 2 ) cos 2 I 

128 \ n V\ n / e2 


+ 15 e 2 (1 - 3 cos 2 I) cos 2g} 


9 / n ®V n ®\ j n 2 sin I (1-3 cos 2 I) {2 (2 + 3 e 2) + 15 e 2 cos 2g} 

e (1 - e 2 ) 2 


32 “ 


9 / n \ t2 u 4 n 2 sin 1(1-3 cos 2 I) 
’ T Kj 22 a 2 (1 - e 2 ) 7/2 


+ — n 2 a 2 sin I cos I { - (2 + 3 e 2 ) + 5 e 2 cos 2g} 


3 . yj 2 
+ — j. R i 


n| sin I cos 1(1-3 cos 2 I 0 ) 


{ 1 — /3 2 (1 + 2 /l - e 2 ) cos 2g } 
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n 2 e cos I (11- 15 cos 2 1 ) 


a (1 - e 2 ) 


sin g 


— J b 4 — — 1 cos 1 { (2 + 3 e 2 ) (3 -7 cos 2 I) - 2 e 2 (4-7 cos 2 I) cos 2g} 
52 a 2 (1 - e 2 ) 7/2 


L5_j s n 2 e cos I {2 (4 + 3 e 2 ) (29- 126 cos 2 I + 105cos 4 I) sin g 
156 a 3 (1 - e 2 ) 972 


- 7 e 2 sin 2 I (7-15 cos 2 I) sin 3g} 


9 / n < t\/ n ®\ 3 n 2 a 2 e cos I 

L 28 \ n * ) V n j /i _ e 2 e 2 


{64-99 e 2 + 3(12- 17 e 2 ) cos 2 I 


+ 15 (2-3 e 2 ) sin 2 I cos 2g} 


i. (Hi) Hi) j b* a! e ^ 1 1 { 2 (7 + 3 e ! ) + 1 S (1 + e 2 ) cos 2g> 

l 6 W/W c(l-e 2 ) 3 


63 / n \ t2 ,4 n 2 e sin 2 I cos I 




2 /1 

a * (1 - e z ) 


— n 2 a 2 e {1-3 cos 2 1-5 sin 2 I cos 2 g} 

8 ® 


3 n£ sin 2 I (1 - 3 cos 2 I ) 

4j,R£ — — fi 2 (2 + 


- e 2 ) cos 2 g 




(1+4 e 2 ) sin g 


15 T , 4 n 2 e 


128 


J 4 b 4 


a 2 (1 - e 2 ) 


2 N 9/2 


{ (4 + 3 e 2 ) (3-30 cos 2 I + 35 cos 4 I) 


- 2 (2 + 5 e 2 ) sin 2 I (1-7 cos 2 I) cos 2g} 


"255 J 5 b 5 " 3 ■ ^ 2 ^ii /2 ^ (4 + 41 e 2 + 18 e 4 )(l-14 cos 2 1 + 21 cos 4 I) sin g 


- 7 e 2 (1 + 2 e 2 ) sin 2 1(1-9 cos 2 I) sin 3 g) 


^-igCfV^Y ol a 2 e 2 /Tr^ sin 2 1 cos j sin 2g 

3g 64 \n*J\ n I e 2 & 


135 / n ®\/ n e 


M/ \ t t 2 e 2 sin 2 I cos I . ^ 

— II— I J 22 b 7 sin 2g 


16 Vn*/\ n 


e (1 - e 2 ) 


2 n 2 


15 

8 




a 2 e 2 sin 2 I sin 2 g 


3 . d2 n l sir > 2 I (1-3 cos 2 I # ) 
+ 4 J 2 R ® i 


/3 2 (1 + 2 / 1 - e 2 ) sin 2 g 



^3 n 2 e sin 1(1-5 cos 2 I) 
a (1 - e 2 ) 5/2 


cos g 
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sin 2g 


15 

32 


J 4 b4 


n 


sin 2 I (1- 7 cos 2 I) 

2 ,i 2 n 7 / 2 

a 2 (1 - e 2 ) 


_ JJL T b 5 — - sin ~— {2 (4 + 3 e 2 ) (1 - 14 cos 2 1 + 21 cos 4 I) cos g 
256 Js a 3 (1 _ e2) 9/2 

- 7 e 2 sin 2 1(1-9 cos 2 I) cos 3g}. 


Now, recalling that 

V" =f f i [ g'l (0* G 'o ] dt 

A 2 g" = / f 2 [g" (t), G''] dt 

a 2 h'' = ; f 3 [ g ; ( t), g;] dt 

A 2 G" = / § [gj (t), G''] dt, 


it is noted that if one performs the integrations, one obtains a small 
divisor (n”) which is of second order. The resulting perturbations will 
therefore be of the first order, and this is not desirable. However, since 
the pericenter moves slowly (about 10 ~ 3 rad/day), these perturbations 
will be of the second order if the time period of integration is of the 
order of hundreds of days. In the following integrals the initial values 
of the Keplerian elements are to be used. For convenience, however, we 
write e for , etc. 
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s f — dt = V— Vlll; 

•'o \ 256 \ n * /\ n / -2 

L \ <w \ / e 


- e 2 cos I [2 + 33 e 2 - (2-17 e 2 ) cos 2 l] 


.unh\ Hi\r b 2 

32 L* U J22 


l2 n 2 sin 2 I cos I /0 . 0 
22 b T ~7 (2 + 3 e ) 


a (1 - e 2 ) e 


117 j n \ t2 U 4 n 2 sin 2 I cos I 1 
~ ~WJ 22 a 3 (1 - e 2 ) 7/2 '8 


+ — J, b 4 ^ 

!28 J4 a 3 (1 - e 2 ) 7/2 


(3-30 cos 2 1 + 35 cos 4 I) (2 + 3 e 2 ) V t 


1 j 3 , , 3 n 2 e sin I (1 - 5 cos 2 I) 
<12 3 a 2 (l-e 2 ) 5/2 


+ i!j. b 5 n 2 e sin I (1 _ 14 cos 2 j + 2 1 cos 4 I) (4+ 3 e 2 )U 

64 Js a 4 (1 - e 2 ) 9/2 J 


(cos g’i - cos g") 


+ 2_( + 945/<y<\ 3 ni e2a 
n" 1 512 l n * j\ n / e 2 


- e 2 cos I sin 2 I 


405 / n <A/ n «\ T i 2 n 2 e 2 sin 2 I cos I 15 2 2 • 2 t 

I — II 1 I„b + — a sin^ I 

128 WA n / a (1 - e 2 ) 2 e 16® 


J 4 b 4 —sin 2 I (1-7 cos 2 I)l(sin2gJ-sin2g") 

1 28 a 3 (1 - e 2 ) ?/2 J 


+ 2|- J b 5 h 2 e 3 sin 3 I (1 _ 9 cos2 I)(cos 3 g » _ cos 3 g") 
!28 a 4 n" (1 - e 2 ) 972 
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I 


dt = ^ (~\(~] ~ a 2 / 1 - e 2 " sin I [2 + 33 e 2 -3(2-17e 2 ) cos 2 I] ' 

31 1 128 U?/W 


-£&)(— \ J„b 2 n2 sin 3- (1-3 cos 2 I) (2 + 3 e 2 ) 


n* U J22 


(1 - e2 ) € 


1 /JL\ j* b 4 n2 Sin I - (1-3 cos 2 I) 
4[n*J J22 a 2 (1 - e 2 ) 7/2 


- — n 2 a 2 sin I cos I (2 + 3e 2 ) 

8 0 

o n 2 sin I cos I 

+ Jh R e (1-3 cos 2 I # ) 


-4ll. b 4 n2skIC ° S i (3-7 cos 2 I) (2 + 3 e 2 ) V t 


32 


a 2 (1 - e 2 ) 


2x7/2 


+ — J — J 3 b 3 — 1 — - cos 1 (11 - 15 COS 2 I) 

n" | 8 J 3 a( i_ e 2)S/2^ 


15 T U5 

+ I28 J3 b 


n 2 e cos I 
a 3 (1 - e 2 ) 9/2 


(29- 126 cos 2 1+105 cos 4 I) (4 + 3e 2 )U 


. (cos g"- cos g") 



e 2 / 1 - e 2 sin 1(1-3 cos 2 I) 



(l-e 2 ) 2 e 


(1-3 cos 2 I) 


1 5 2 2 2 T t 

+ — n a * e* sin I cos I 
16 ® 


126 



15 - , 4 n 2 e 2 sin I cos I 
32 4 a 2 (1 - e 2 ) 7/2 


(4-7 cos 2 I) 


3 n 2 sin I cos I (1 - 3 cos 2 I # ) (3 2 (1 + 2 /l - e 2 ) 

7 


• ( si n 2 g" - sin 2 g") 


- 35 y b 5 


n 2 e 3 cos I 


256 5 n" a 3 (1 - e 2 ) 9/2 


sin 2 I (7- 15 cos 2 I) (cos 3 g'j- cos 3 gJJ) 


*^idt = /^V— V 

|128 \n*/\ n / (i_ e 2^ 1 / 2 e 2 


a z e cos I 


[64-99 e 2 + 3 (12-17 e 2 ) cos 2 I] 


9 / ri^\ /n^\ 9 « 

( Z \ I 1 K2 n e c : «2 T t /7 . 3 ~2> 

8 W\ n 7 22 (1 - e 2 ) 3 e 


sin 2 I cos I (7 + 3 e 2 ) 


63 / n \ T 2 k 4 n 2 e 




4 WJ 22 a 2 (1 - e 2 ) 972 


sin 2 I cos I — — n 2 a 2 e (1 - 3 cos 2 I) 


1 5 , , 4 n 2 e 
128 4 a 2 (l-e 2 ) 9/2 


(3-30 cos 2 11+35 cos 4 I)(4+3e 2 ) L t 


+ J, b 3 ° 2 Sln l — (1-5 cos 2 I) (1 + 4 e 2 ) 


n" ] 8 J 3 „ ,1 _^V2 


a (1 - e 2 ) 


+ J s b 5 — — S * n ^ 1/2 (1-14 cos 2 1 + 21 cos 4 I)(4 + 41 e 2 + 18 e 4 ) 

a 3 (1 - e 2 ) 


.(cos g"- cos g") 


127 


d / n a~ e 

5 \ n *y \ n / e 2 (1 - e 2 ) 1/2 


cos I sin 2 I (2 - 3 e 2 ) 


-^V— ) J„b 2 h 2 *— sin 2 I cos I(l + e 2 )+i|-n| a 2 esin 2 I 




R 2 — 5 — sin 2 I (1-3 cos 2 I 0 )/3 2 (2 + /l _ e 2 


J 4 b 4 n 6 — sin 2 I (1-7 cos 2 I)(2 + 5 e 2 ) Wsin 2 g",- sin 2g"; 

a 2 (l-e 2 ) 9/2 J 


n 2 e 2 


5 u ~ 2N 11 / 2 

n a 3 (1 - e 3 ) 


sin 3 1(1-9 cos 2 I) (1 + 2 e 2 )(eos 3g" - co 


3 _ .a n 2 e sin I ,, i 15 T L s n 2 e sin I 

8 Jsb 5cos I)+ra,s 


(1-14 cos 2 I + 21 cos 4 I) (4 + 3 e 2 )> (sin g'{- sin g^) 


f 135 /M/^V — 
^ + 128\^n*J\nJ e 2 


a 2 e 2 


? 2 sin 2 I cos I 


RV— 1 J 22 b 2 ° 2 e2 r sin 2 I COS I 


32 \<\ n ) 22 ^(l-e 2 ) 2 


+ — n 2 a 2 e 2 sin 2 I 


_i; r 2 sin 2 I (1 - 3 cos 2 I ffi ) /3 2 (1 + 2 /l - e 2 



1 5 T , 4 n 2 e 2 


+ ^J 4 b 4 


-sin 2 I (1 - 7 cos 2 I) l(cos 2 g"- 


64 ‘' 4 ~ a 2 (1 - e 2 ) 772 ~ *' jv— 61 -cos2g") 


35 T , c n 2 e 3 


+ h b 5 


256 5 n" a 3 (1 - e 2 ) 9/2 


sin 3 1(1-9 cos 2 I) (sin 3gj - sin 3 g") 


34. Summary of the Development 

The short-period terms are given by 


/ = /' - 


BS, 


dV 


, 3S 2 

g = g 

BG 


h = h' - 


BS, 

BH' 


L = L' + 


B^ 

B/ 


G = G' + 


bs 2 

3 g 


H = H' + 


BS, 

Bh 
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Long-period terms are given by 


. Hi 

. Hi 

_ Hi 

BL" 

BL" 

BL" 

bs; 

_ Hi 

bs; 

1 ^ 

1 ro 

BG" 

BG" 

_ Hi 

BS* 

. Hi 

BH" 

BH 7 

BH" 


G' 


G" 


Hi 

Bg* 


Hi + HI 

Bg' + Bg' 


, BS BSl Bs; 

H' = H" + + — + - 

Bh' Bh' Bh' 


Additional terms are obtained from 


J (-- - S) - 

fBF" 

A. l" = - — dt 

2 J BL" 


r -ii = Aj r + a 2 v 
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g" - go = A, g" + A 2 g- 

A - h " = - 1 ("« + P) dt 

r bf 3 

" = - — dt 

J BH" 


a 2 h 


h" - h" = Aj h" + A 2 h" 


L" - L" = 0 


Aj G" = 


r BF_" 


a 2 G 


„ _ f bf; 

J Bg 


dt 


dt 


G" - G" = Aj G" + A 2 G" 


H" - H" = 0 . 


Then we have 


AL = L - L' 


AL' = 0 


AL" =0 


A G" = G" - G" 


A G = G _ G' A G' = G' - G" 

AH=H-H' AH' = H' -H", and AH" rO. 


Note that in the expressions for AjG, A -t , g, and A h the variables 
G" + A 2 G" and gg + A 2 g" should be used in computing v 2 , <9" 2 ,C, etc. 

Thus, the perturbations in the Keplerian elements a, e, and I are 
obtained by 


a 


+ 


2a' 


AL 

V 
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e 


= e 


+ 


1 - e 


' 2 


( 


al 

V 



Also , 


I 


I' + cot I' 




a 


n 


e 


# 


= e 


1 - e" 2 AG' 
e" G" 


1 - e" 2 /AG'\ 2 
2 e" 3 \ G" / 


I' 


/AG' AH'X cos I* /AG'\ /AH'\ 

ycF ~ ir ) 4 sin 3 1" yG* - ) yiP ” ) 


1 cos T " 
2 sin 3 I" 


(1 + sin 2 I") 



1 cos 3 I" /AH'Y 
“ 2 sin 3 I" \H" ) ’ 

where in the squares and products of AG' and AH' terms of the 3rd 
order or superior should be neglected. Finally, 


a 


N 


a 


u 

0 


« 


e 




I" = I" + cot I 


n 

0 


AG" 


35 . Position and Velocity 


e/0; 1/0°, 180° 

From the elements L, G, H, / , g> h, one can obtain the 
coordinates and the components of velocity as follows: 

Obtain a , e , and I from the following equations: 
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a = L 2 /m 


I = arc cos - C0° < I < 180°) 

G 


Now solve Kepler's equation 


E - e sin E r / 


to obtain E. 


Then compute r from 


r = a (1 - e cos E) 


Next obtain f from 


cos f = — (cos E - e) 


aG . r 
sm f =y — Sln E. 
^ L 


Now compute 


A = a [cos g cos h - sin g sin h cos I] 


B = a — [- sin g cos h - cos g sin h cos I] 

x L 


A y =: a [cos g sin h + sin g cos h cos I] 


B y = a— [- sin g sin h + cos g cos h cos I] 


A =a sin g sinl 


B z = a— cos g sin I. 
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Then, 



Since the system is rotating, 
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37. Appendix A 


The elliptic integrals that are needed for this section may be found 
in Reference 3. They are 


f dx 

Jxj (x - V 2 ) 


2 ) /Q(x) 


2 . 


r x 

J xl x 


dx 


X[ X /Q(x) 


f 42- 

JxJ fl - x) V 


(1 - x) v'Q(x) 


!I 


dx 


i x 2 V/QOO 


5. 


£ 


dx 


x' x 3 v Q(x) 


6. 


7. 


8 . 


9. 


10 . 



dx 

(x - V 2 ) / xQ(x) 
dx 

(1 - x) / xQ(x) 


dx 

/ xQ(x) 


dx 

x / xQ(x) 
dx 

x 2 /xQ(x) 


where Q(x) = (x - xj) (x - x‘) (x - x' 3 ), and x' > x' - x > x[ > 0. 


1-5 


Let 


x - x, 

2 - 1 


sin 2 4 1 = 


I I 

x 2 - x r 


k 2 _ X 2 - X 1 
X 3 - X 1 
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We then have 


r 

J ^ i 


dx 


i (p - *) m ^QOO (p 


-h— r 

- xj) m /xj - xj 0 


du 

(1 - a 2 sn 2 u)' 


where 


x' - xj 

snu = sin i p, a 2 = 

p - x! 


For m = 1: 


i 


dx 


xj (p - x) v/Q(x) 


2 I du 

( P -x;)v^T^[ -l 1 - a2 s " 2u 


n(0, a- 2 , k). 


(p - x{) vx 3 - x' 


Thus, 


I 


dx 


II ('Z', a 2 , k), a : 


x l (x - V 2 ) /Q 00 ( X ; _ v 2 ) /x* - x ; 


2_ X 2- X I 

l , 2 - xj 


I 


dx 


n(0, a 2 , k) , 


x # — x # 
2 X 2 X 1 
a 46 = - 


[ X /Q(x) x | /x^ - x| 


I 


dx 


xj (1 - x) v Q(x) - 7~7== 

(1 - Xj) /x 3 - Xj 


n a 2 , k) , 


2 X 2 " X I 

a = 


i- x i 
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For m = 2: 


f dx 2 f U du 

•l; x^/qw -l a-* 2 sn 2 u ) 2 


t I 

x 2 - Xj 


X l 2 17 X 3 - X 1 


2 (a 2 - 1) (k 2 - a 2 ) 


a 2 E(u) + (k 2 - a 2 ) u + (2 a 2 k 2 + 2a 2 - a 4 - 3k 2 ) 11(0, a 2 , k) 


a 4 sn u cn u dn u 
1 - a 2 sn 2 u 


where u = F(i p, k ) 
and E(u) = E (\p, k) 
For: m = 3: 


f dx _ 2 f du 

X 3 /Q00 x ;3 i - “ 2 s " 2 u ) 3 


x' - x' 
x 2 X 1 


x l 3 ' /x 3- x l 


where 


— 2T7T2 2l k2 u + 2 ( a2 k 2 + a 2 - 3k 2 ) nOA, a 2 , k) 

4(1 - a 2 ) (k 2 - a 2 ) 


+ 3 (a 4 - 2a 2 k 2 - 2a 2 + 3k 2 ) V + 


a 4 sn u cn u dn u 


(1 - a 2 sn 2 u) 2 
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6 & 7 


Let 


sin 2 \p = 




x'(x» - xp 
x '(x' - X')’ 


We have 



* 1 


dx 

(p - x) m v'xQ(x) 


2 

(P~ xj) m /x'(x^ - x') 



(1 - a 2 sn 2 u) m 
(1 - a 2 sn 2 u) m 


where 


snu = si n \p t 


P(x^ - x|) 

x^(p - X' ) ' 


For m = 1, p^O: 


i: 


dx 


xj (P “ X) VxQ(x 


_ r u 

2 f 1 - a 

) (p - xp a/x'(x' - xp Jo 1 - a : 


2 2 
sn z u 


2 2 
sn" u 


du 


p/x'(x' - xj) 


u + r n(0, a* k) 

P - x i 


Thus, 


J .x 

dx 

xj (x - ^ 2 )/xQ(x) 


- 2 



u 



n<0, a 2 , k) , 


v 2 (x* - xj) 
x'(v 2 - xj) 
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i: 


dx 


(1 - x)/xQ(x) /x'(x' - *;") 


u + 


1 - 


n (0. k) 


a 3 = 


- X, 


x'(! - xj) 


8-10 


Again, let 


sin 2 0 = 



We have 



x 3( X 2 * X i> 

X 2< X 3 - x l>’ 


a 2 = 



f X dx 
Jxj x m / xQ(x) 




/xlc 


_r 

4 - xj) J 0 


(1 - a 2 sn 2 u) ro du. 


For m = 0: 


J mX 

dx 

xj / xQ(x) 


2 

/ x 2( x 3 - X D 


For m = 1 : 


f* dx 2 f U „ 2 2 n 

— = — I (1 - or sir' u) 

J x j x/xQ(x) xj /x'(x' - x j ) Jo 


du 


x ;/x'( x ; - x ;> 


[(k 2 - a 2 )u + a 2 E(u)] 


X 1 X 3 A,(x' - xj) 


— [xj u + (x' - x{)E(u)]. 
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For m = 2: 


r. 


dx 


x 2 / xQ(x) Xj 2 / Xj (x' - xj) 


f 


(1 - a 2 sn 2 u) 2 du 


\ — — [(3k 4 - 6a 2 k 2 + 2a 4 +k 2 a 4 )u 

xj 2 / x' (x' - x[) 3k 4 


+ 2 (3 a 2 k 2 - a 4 - k 2 a 4 ) E(u) + a 4 k 2 sn u cn u dn u]« 
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3 8. Appendix B 


The integrals that are needed for this section may be found in 
Reference 8. They are 




1 . 



The second of these integrals may be found in Part 1 . 
3 . 
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4 . 


(x - v 2 )(x - l)/x(x - X-) V 2 - 1 [J(x - V 2 )/x( 


If d > + f — d x . . 1 

]J(x - v 2 )/x(x - x 3 ) J(1 - x)v^c(x - X 3 )^ 


2v 2 (v 2 -x 3 ) + (2v 2 -x 3 )(x - v 2 )~2 Sx(x-x 3 )(y 2 -x.)v 


( v 2 > x 3 ) 


(x - v 2 ) / x(x-x 3 ) ^x 3 -v 2 )v 2 


-2v 2 (x 3 - z/ 2 ) + (2v 2 - x 3 )(x - v 2 ) 
2 /x(x - x, )(x, — v 2 ) Z 2 


(v 2 < x,) 


(1 - x)/x(x - X 3 ) /I - X 3 


1 2(1 -X 3 ) + (x 3 - 2)(1 -x) + 2 /x(x -x 3 )(l - X 3 ) 

log 


r dx 
J x(x - 1) /? 


x(x - x 3 ) 


r dx_ 

J x /x(x - 


x 3 ) J (1 - x) /x(x - x 3 ) 


f— . 

J x/x(x - X-) 


^x(x - X,) 


The second of these integrals may be found in Part 4. 


(x - l)/x(x - x-) /I - x 


1 2(1 - x 3 ) + (x 3 - 2)(1 - x) + 2 /x(x -x 3 )(l - x 3 ) 

log 


1 - X 


j] dx _ |* dx dx + j* 

Jx 2 (x - 1 )/x - x 3 Jx 2 /x - x 3 J x/x - x 3 J (x - 


l)i/?c - x, 


x vx - x 3 x 3 x 2 


1_ f dx 
X, » Xi/x - X, 
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The integrals 


I 


dx 


and 


(x - 1) /x - x^ 




dx 


/x - X, 


may be found in Parts 1 and 2, respectively. 


8 . 


1 

II 

X 

i 

dx 

dx 

* dx 

+ 

- dx 

J x 3 (x - l)/x - X 3 

1 ! 
) x J /x - x 3 v 

| ? / - " ' 
ix> x - x, 

' 3 * 

! v/x - x, J 

(x - l)/x - X 3 


I 


dx 

X 3 / X - X 3 


yx=ir 3 3 
2x 3 x 2 4x 3 



dx 


The integrals 


f d * , 1 

dx | 

and 

dx 

J x 2 /x - X 3 

|(x - 1) i/x - x 3 

1 x/x - x 3 


may be found in Parts 7,1, and 2, respectively. 


9. 


! 


dx 

(1 - x) 2 / x(x - x 3 ) 


/x(x-x 3 ) t 2 -», (• dx 

(1 - x 3 )(l - x) 2(1 - x 3 ) J (1 - x) / x(x - X 3 ) 
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